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ABSTRACT 

We use imaging from the Next Generation Virgo cluster Survey (NGVS) to present a comparative 
study of ultra-compact dwarf (UGD) galaxies associated with three prominent Virgo sub-clusters: 
those centered on the massive, red-sequence galaxies M87, M49 and M60. We show how UGDs can 
be selected with high completeness using a combination of half-light radius and location in color-color 
diagrams {u'^iKs or u^'gz). Although the central galaxies in each of these sub-clusters have nearly 
identical luminosities and stellar masses, we find large differences in the sizes of their UGD populations, 
with M87 containing ^ 3.5 and 7.8 times more UGDs than M49 and M60, respectively. The relative 
abundance of UGDs in the three regions scales in proportion to sub-cluster mass, as traced by X-ray 
gas mass, total gravitating mass, number of globular clusters, and number of nearby galaxies. We 
find that the UGDs are predominantly blue in color, with ^ 85% of the UGDs having colors similar 
to blue GGs and stellar nuclei of dwarf galaxies. We present evidence that UGDs surrounding M87 
and M49 may follow a morphological sequence ordered by the prominence of their outer, low surface 
brightness envelope, ultimately merging with the sequence of nucleated low-mass galaxies, and that 
envelope prominence correlates with distance from either galaxy. Our analysis provides evidence that 
tidal stripping of nucleated galaxies is an important process in the formation of UGDs. 

Subject headings: galaxies: dwarf; galaxies: individual (M87, M49, M60); galaxies: nuclei; galaxies: 
star clusters: general; galaxies: structure; globular clusters: general 
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1. INTRODUCTION 

For many decades, it was widely believed that galaxies 
and star clusters were completely unrelated populations 
in the manifold of stellar systems, characterized by dis¬ 
tinct stellar populations, structural and dynamical prop¬ 
erties, and, presumably, formation histories (e.g., Kor- 
mendy 1985; Burstein et al. 1997; see also the overview 
of Willman & Strader 2012 and references therein). How¬ 
ever, in the late 1990s, multiple spectroscopic surveys of 
the Fornax cluster uncovered an apparently new type 
of compact stellar system with properties that seem¬ 
ingly bridged the gap between compact, low-mass galax¬ 
ies and globular clusters (GGs). First, Hilker et al. 
(1999) used the 2.5 m Dupont telescope to carry out 
a spectroscopic survey of stellar systems in the vicin¬ 
ity of NGG1399 and reported the discovery of two For¬ 
nax members so compact that they were barely resolved 
in ground-based images. Three other compact cluster 
members were identified by Drinkwater et al. (2000) and 
Phillipps et al. (2001) in an AAT/2dF Fornax spectro¬ 
scopic survey. With absolute magnitudes in the range 
— 13.8 < Mb ^ —11.6, these objects have luminosities 
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comparable to those of the faintest dwarf galaxies then 
known in the Local Group, but effective radii about an 
order of magnitude smaller. As a result, these rare sys¬ 
tems came to be known as “ultracompact dwarf’ galaxies 
(UCDs). 

Soon afterwards, UCDs were discovered in the Virgo 
cluster using a combination of HST imaging and Keck 
spectroscopy (Hasegan et al. 2005). Additional Virgo 
UCDs were later identified in an AAT/2dF spectroscopic 
survey (Jones et al. 2006). More recent studies have in¬ 
creased the number of known UCDs in Virgo, though 
mainly in the immediate vicinity of M87 (e.g., Brodie 
et al. 2011). To date, UCDs have been found in a wide 
range of environments, including the Centaurus (Mieske 
et al. 2007), Hydra (Wehner & Harris 2007; Misgeld 
et al. 2011), Coma (Price et al. 2009; Madrid et al. 2010; 
Chiboucas et al. 2011), Antlia (Caso et al. 2013) and 
Perseus (Penny et al. 2014) clusters, as well as in prox¬ 
imity to some group and field galaxies (Evstigneeva et al. 
2007; Madrid 2011; Hau et al. 2009; Norris & Kannappan 
2011). 

Generally speaking, the UCD samples amassed by 
these surveys have continued to close the once-prominent 
gap between GCs and “normal”, low-mass galaxies, fur¬ 
ther blurring the distinction between these families of 
stellar systems.Still, the role of selection effects 
in shaping our perceptions of “families”, “gaps”, and 
“trends” among stellar systems is a matter of concern. 
Existing UCD catalogues are heterogeneous in nature, 
having been assembled from a patchwork of surveys with 
different biases and completeness limits. This is under¬ 
standable because UCDs are notoriously difficult to find 
and study: i.e., most are too small to be resolved in typi¬ 
cal ground-based images, rendering spectroscopic surveys 
the obvious route forward (Hilker et al. 1999; Drinkwa- 
ter et al. 2000; Jones et al. 2006; Mieske et al. 2007; 
Evstigneeva et al. 2007; Adami et al. 2009; Gregg et al. 
2009). Unfortunately, spectroscopy is both time consum¬ 
ing and inefficient: for instance, the 2dE survey of Jones 
et al. (2006) yielded just nine UCDs from a radial veloc¬ 
ity sample of 1501 objects, for an overall success rate of 
- 0 . 6 %. 

In principle, high-resolution Hubble Space Telescope 
(HST) imaging can circumvent this problem by allowing 
UCDs in the local universe to be resolved directly and 
separated from GCs on the basis of size and luminosity 
(Mieske et al. 2004b; Hasegan et al. 2005; Blakeslee & 
Barber DeGraaff 2008; Hau et al. 2009; Madrid et al. 
2010; Madrid 2011; Strader et al. 2012; Caso et al. 
2013; Norris et al. 2014). Although this approach can 
improve discovery efficiency, candidates must still be 
confirmed spectroscopically because contamination from 
background galaxies can be significant, particularly at 
faint magnitudes. The foremost concern, though, is 
HST’s small field of view, which severely limits sky cover¬ 
age and hence spatial completeness in the nearest groups 
or clusters: e.g., the search for UCDs in the Virgo cluster 
by Hasegan et al. (2005), undertaken as part of the ACS 
Virgo Cluster Survey (ACSVCS; Cote et al. 2004), cov- 

A working definition of UCDs has recently been proposed by 
Willman & Strader (2012), who take UCDs to be systems with 
sizes 10 < r/i < 100 pc and absolute magnitudes —13 < My < 
—9. However, a variety of such definitions exist, and it is worth 
remembering that all such definitions remain subjective in nature. 


ered one hundred sightlines but just ^0.3% of the total 
cluster area. 

More than 15 years after their discovery, the origin of 
UCDs remains surprisingly obscure. A number of scenar¬ 
ios — with varying degrees of theoretical/numerical un¬ 
derpinning — have been proposed: UCDs could be: (1) 
otherwise normal objects within the high-luminosity tail 
of the GC luminosity function; (Mieske et al. 2002, 2012); 

(2) end-products of the aggregation of young massive star 
clusters formed during the interaction of gas-rich galaxies 
(Eellhauer & Kroupa 2002, 2005; Briins & Kroupa 2012); 

(3) the remnants of nucleated dwarf galaxies that have 
been stripped, or “threshed”, during tidal interactions 
(Bekki et al. 2001; Drinkwater et al. 2003; Paudel et al. 
2010; Pfeffer & Baumgardt 2013; Seth et al. 2014; Janz 
et al. 2015). 

Of course, these scenarios need not be mutually exclu¬ 
sive, and it is entirely possible that multiple formation 
pathways exist, as noted by Hasegan et al. (2005) and 
others (Mieske et al. 2006a; Brodie et al. 2011; Chilin- 
garian et al. 2011; Da Rocha et al. 2011; Norris & Kan¬ 
nappan 2011). Eor instance, Mieske et al. (2012) ex¬ 
amined the specific frequency of UCDs in different en¬ 
vironments and showed that they seem to match those 
of GCs very well, suggestive of a direct GC-UCD link. 
Other studies found some UCDs to be surrounded by 
faint halos (Hasegan et al. 2005; Blakeslee & Barber 
DeGraaff 2008) or even asymmetric extensions (Richtler 
et al. 2005; Brodie et al. 2011) that could be explained by 
nucleated dwarf progenitors. There is also evidence that 
mergers of super star clusters has played a role in the 
formation of some young massive objects that could be 
analogs for some UCD progenitors (e.g., Maraston et al. 
2004). Thus, the immediate tasks are to decide which 
mechanisms have been involved in UCD formation, and 
to assess their relative importance as a function of time 
and environment. 

A large and carefully selected sample of UCDs — 
preferably identified in a survey with a well character¬ 
ized selection function — is needed to better understand 
the origin of these puzzling systems. As the nearest rich 
cluster of galaxies (Mei et al. 2005; Blakeslee et al. 2009), 
Virgo is an obvious choice for a comprehensive UCD sur¬ 
vey. At a distance of ^16.5 Mpc, a UCD in Virgo with 
a typical half-light radius of ^ 20 pc (corresponding to 
^ 0725) can, as we will show below, be resolved in high- 
quality ground-based images. Moreover, the cluster itself 
is quite rich, containing ^ 2000 cataloged member galax¬ 
ies of virtually all morphological types (Binggeli et al. 
1987), tens of thousands of GCs, and large numbers of 
nucleated galaxies (Binggeli et al. 1987; Cote et al. 2006; 
Jordan et al. 2005; Harris 2009a), making it possible to 
compare the properties of these stellar systems within a 
single, homogenous survey. 

The Next Generation Virgo cluster Survey (NGVS) is a 
large program carried out with the 3.6 m Canada-Erance- 
Hawaii Telescope (CEHT) between 2008 and 2013. The 
project has been described in detail in Eerrarese et al. 
2012 (Paper I). Briefly, the survey used the MegaCam 
instrument to perform panoramic imaging of the Virgo 
cluster, from its center to virial radius, in the u'^giz fil¬ 
ters (^ 100 deg^). Eor a subset of the survey area, r- 
and ATg-band imaging is also available (see, e.g., Munoz 
et al. 2014). Especially important for the study of UCDs, 
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the NGVS image quality is uniformly excellent, with a 
median i-band seeing of FWHM 0'/54 (the first and 
third quartiles are ^ 0'/52 and 0'/56, see also Figure 8 
in Ferrarese et ah 2012). Overall, the depth, areal cover¬ 
age, and image quality of the survey, as well as the avail¬ 
ability of extensive color information, make it an ideal 
resource for the study of compact stellar systems like 
UCDs. We will identify UCD candidates using NGVS 
multi-wavebands data and show that this photometry- 
based sample is very clean. Because the survey spans a 
wide range of local density, its also provides an opportu¬ 
nity to investigate the role played by environment in the 
formation and evolution of UCDs. 

In this paper, we use NGVS imaging to carry out a sys¬ 
tematic study of UCDs in the Virgo cluster. We focus on 
three important, high-density regions within the cluster: 
the sub-clusters centered on the massive, red-sequence 
galaxies M87, M49 and M60. In future papers, we shall 
extend our analysis to include the entire cluster, includ¬ 
ing regions of low density. Other papers in the NGVS 
series related to the topics considered here include stud¬ 
ies of the cluster-wide GC populations in Virgo (Durrell 
et al. 2014), the abundance and dynamical properties of 
star clusters, UCDs and galaxies in the cluster core (Zhu 
et al. 2014; Zhang et al. 2015; Grossauer et al. 2015), 
the structure and dynamics of compact galaxies in Virgo 
(Guerou et al. 2015) and the physical classification of 
stellar and galactic sources using optical and NIR imag¬ 
ing (Munoz et al. 2014; Raichoor et al. 2014). Zhang 
et al. (2015) is especially relevant to this paper, in that 
it presents a detailed kinematic analysis of UCDs and 
GCs in the core of the cluster. 

This paper is structured as follows. We give a brief de¬ 
scription of the NGVS imaging and data reduction pro¬ 
cedures in §2.1 and §2.2, while a discussion of our UCD 
selection methods is given in the remainder of §2. In §3, 
we examine and compare the properties of UCDs in the 
M87, M49 and MOO regions, paying particularly close at¬ 
tention to the M87 region where we have the benefit of 
both TCg-band imaging and extensive spectroscopic cov¬ 
erage. The implications of our findings are given in §4, 
while §5 presents our conclusions and outlines some pos¬ 
sible directions for future work. 

2. OBSERVATIONS 

2.1. Survey Overview and Data Reduetion 

The NGVS is a deep optical imaging survey in the 
i and z bands, with additional partial coverage in r, 
carried out with MegaCam, the 340 megapixel camera 
on CFHT (Boulade et al. 2003).^^ The survey covers an 
area of U = 104 deg^ inscribed within the virial radii of 
the Virgo A (R 200 = 1-55 Mpc) and Virgo B (R 200 = 0.96 
Mpc) subclusters (Ferrarese et al. 2012), which are cen¬ 
tered on M87 and M49, respectively. Beginning with 
the NGVS “pilot program” (see Ferrarese et al. 2012), 
images were collected during six consecutive observing 
seasons, 2008 through 2013. Each MegaCam pointing 
(i.e., NGVS field) measures roughly 1° x 1° on the sky. 
There is some overlap between fields, so that the 104 deg^ 
survey area is covered by 117 distinct NGVS pointings. 
Each NGVS field is assigned a pair of ordered numbers 

Outside of the M87 region, r-band imaging is available at only 
a fraction of the full survey depth, ranging from 1374 to 4461 sec. 


according to its location within the survey grid. Eield 
(+0,+0) is in the cluster center and includes M87; the 
field numbers grow with increasing right ascension and 
declination. In addition, four background fields were also 
observed (see Table 1 of Eerrarese et al. 2012). The back¬ 
ground fields are located well beyond the survey bound¬ 
aries (offset from M87 by 16°), at about three virial 
radii from the center of the A subcluster and at Galactic 
latitudes corresponding to the highest and lowest values 
spanned by the main survey. 

Exposure times in seconds for the five filters are: 

Tu* = 11 X 582 = 6402 

Tg = 5 X 634 = 3170 

Tr = 7 x 687 = 4809 (1) 

Ti = 5 x 411 = 2055 

= 8 X 550 = 4400 

Different dithering patterns were used for the different fil¬ 
ters with the exception of the g and i bands which shared 
an identical strategy. In the ^-band, the limiting magni¬ 
tude for point sources is 25.9 mag (lOa) while the surface 
brightness limit for extended sources is jig ^ 29.0 mag 
arcsec“^ at 2a above the mean sky level. This means 
that the NGVS is roughly three magnitudes deeper than 
either the Sloan Digital Sky Survey (SDSS; York et al. 
2000) or the Virgo Gluster Gatalog (VGG; Binggeli et al. 
1985). In addition, the NGVS image quality is excellent: 
the EWHM for every image, in all filters, is always bet¬ 
ter than 1". The best seeing conditions were reserved 
for the i band, where the EWHM never exceeded 076, 
with a median of 0754. The g band images have longer 
exposure time and the larger EWHM with a median of 
0780 (Eerrarese et al. 2012). 

The pre-processing for all NGVS images includes bad 
pixel masking, overscan and bias subtraction, fiat field¬ 
ing, fringing correction (which is necessary only for the 
i and z bands) and scattered light correction. The in¬ 
dividual images are then corrected for a global back¬ 
ground map and stacked using the MegaPipe pipeline 
(Gwyn 2008) with the adoption of an artificial skepti¬ 
cism method (see Eerrarese et al. 2012 for details). The 
final stacks are then photometrically and astrometrically 
calibrated to SDSS images and scaled to have photomet¬ 
ric zero point of rriAB = 30. 

Eor the Virgo core region centered on M87 (a 3.62 deg^ 
survey that is comprised of the NGVS+0+0, NGVS+0+1, 
NGVS-1+0 and NGVS-1+1 fields), a deep Kg band imag¬ 
ing survey (Next Generation Virgo cluster Survey - In- 
fraRed, hereafter NGVS-IR) was also carried out with 
GEHT, using the WIRGam instrument (see Munoz et al. 
2014 for more details). The NGVS-IR has a limiting 
magnitude of Kg ^ 24.4 mag (5cr) for point sources and 
a surface brightness limit of fiKs — 24.4 mag arcsec“^. 
The seeing was always better than 077, with the median 
seeing of the stacked mosaic being 0754. The NGVS-IR 
raw images were processed using the ‘I‘iwi processing 
pipeline v2.0 with cosmic rays removed and satu¬ 
rated pixels corrected. After the sky removal, the im¬ 
ages were stacked together using the SWARP software 
package (Bertin et al. 2002). Eull details on NGVS-IR 
reductions can be found in Munoz et al. (2014). 

see http://www.cfht.hawaii.edu/Instruments/ 
Imaging/WIRCam/IiwiVersion2Doc.html 
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2.2. Catalog Generation and SExtraetor Photometry 

To construct an NGVS compact-source catalog, 
SExtraetor (Bertin & Arnouts 1996) was first run on 
each NGVS field in double image mode. Objects were 
detected in the ^-band images using SExtraetor pa¬ 
rameters: DETECT_THRESH=1.5, ANALYSIS_THRESH=5.0, 
DETECT_MINAREA=2.0, DEBLEND_NTHRESH=32.0 and 
DEBLEND_MINC0NT=0.002, and then measured in each of 
the E'griz images. More than 25 million objects were 
detected and measured across the 117 NGVS fields. 
All the objects were detected on the MegaPipe global 
background stacked images, which have the highest 
photometric accuracy for point-like sources such as stars 
and GGs. Note, however, that this catalogue is not 
suitable for investigations of extended objects, especially 
low surface brightness galaxies. Euture papers in this 
series will discuss the identification and characterization 
of extended sources in the survey. 

By contrast, UGDs at the distance of Virgo are both 
bright and compact (see, e.g., Hasegan et al. 2005 for an 
early demonstration of this based on HST/AGS imaging) 
so we expect to detect virtually every UGD in the cluster. 
As discussed below, the main obstacle facing a systematic 
study of UGDs in Virgo is the high level of contamination 
from stars, compact background galaxies and spurious 
features such as asterisms or blends. 

2.3. Measurement of Half-Light Radii: Methodology, 

Simulations and Consisteney Cheeks 

While UGDs are often taken to have half-light radii, 
rh, larger than 10 pc, the exact dividing point with GGs 
is somewhat arbitrary. Nevertheless, the great majority 
of GGs are smaller than rh ^ 10 pc; in the Virgo cluster, 
which lies at a distance of 16.5 Mpe (Mei et al. 2005; 
Blakeslee et al. 2009), the average GG half-light radius is 
2.7 ±0.35 pc ^ 07035 (Jordan et al. 2005, similar results 
can be found in Webb et al. 2013). By contrast, the 
half-light radii of known UGDs span a range, from 11 
pc to 93 pc with a median of ^ 20 pc or 0725 (Hasegan 
et al. 2005; Evstigneeva et al. 2008). Virgo galaxies are 
much larger than GGs or UGDs: i.e., the effective radii 
of even the most compact galaxies in Virgo are larger 
than 150 pc ^ 2" (Eerrarese et al. 2006). Thus, a typical 
GG in Virgo can only be resolved at HST resolution, but 
a typical UGD can be resolved in high-quality ground- 
based images, like those from the NGVS. Note that the 
MegaGam pixel size is 07187, which corresponds to ^ 
15 pc at this distance. 

In order to measure structural parameters — includ¬ 
ing half-light radii and total magnitudes — for the many 
thousands of GG candidates in the AGSVGS, Jordan 
et al. (2005) developed the KINGPHOT software pack¬ 
age. Briefly, KINGPHOT fits — inside a fitting radius, 
rfit — PSE-convolved King (1966) models to each GG 
candidate and then determines structural parameters via 
minimization. Needless to say, reliable PSE determi¬ 
nation is essential for the successful application of this 
code. In the case of the AGSVGS, the KINGPHOT mea¬ 
surement limit on r^ was found to be ^ 1 pc, which 
is about l/8th of the instrumental PSE for AGS/WEG 
(EWHM ^ 071). The uncertainty on rh was also found 
to grow when rh became comparable to rfif 

In the present study, KINGPHOT has been adapted to 


measure the half-light radii for UGDs and bright GGs in 
the NGVS images. Not surprisingly, the main challenge 
in an application of KINGPHOT to ground-based images 
lies in the estimation of the PSE. Although the NGVS 
image quality is superb for a ground-based survey, the 
PSE is nevertheless much broader than that of HST (by 
about a factor of five in the i band). In addition, ground- 
based PSEs can change from field to field, and may vary 
significantly across an individual MegaGam pointing. 

Gomplete details on the construction of PSEs in the 
NGVS will be presented in Gwyn et al. (2015). Briefly, 
candidate PSE stars were first identified using SExtraetor 
and DAOphot (only stars identified by both SExtraetor 
and DAOphot were used). Additionally, stars within a 
distance of 40 pixels from other bright objects were dis¬ 
carded. The PSE was then generated using the DAOphot 
routine psf with a bivariate Gaussian and allowed to vary 
across the field with second order variations. In each 
NGVS field, the residual images for the PSE stars were 
checked to ensure there were no systematic trends with 
magnitude or position, and KINGPHOT was run only 
on bright point sources with high signal-to-noise (S/N) 
ratios. Measurements in the g and i bands (which have 
the highest S/N and the best image quality) were made 
independently. Our final rh values represent weighted 
averages of the measurements in the g and i bands. 

Eor the catalog described in §2.2, a set of aperture 
magnitudes (within diameters of 3, 4, 5, 6, 7, 8, 16 and 
32 pixels)were measured for all objects. Because the 
PSEs differ from field to field, the magnitude within a 
given aperture does not represent the same percentage of 
total flux. Accordingly, we corrected the aperture mag¬ 
nitudes to infinite aperture in a two-stage process. Eirst, 
a sample of bright (but unsaturated) stars were selected 
in each field and their 16-pixel aperture magnitudes cor¬ 
rected to infinity (calibrated to SDSS PSE magnitudes). 
Then, the other aperture magnitudes (3, 4, 5, 6, 7 and 8 
pixels) were corrected to a 16-pixel aperture. 

After the application of these corrections, we have a ho¬ 
mogeneous database of point-source photometry across 
the entire NGVS survey area, all calibrated against the 
SDSS. This catalog also accounts for the small but mea¬ 
surable spatial variations in the PSE within individual 
fields, thereby tightening the stellar color-color and color- 
magnitude diagrams and aiding in the selection of targets 
(see, e.g., Eigure 1 of Durrell et al. 2014). Eigure 1 shows 
the aperture correction map made by stacking individual 
MegaGam pointings. The color coding highlights regions 
of negative and positive deviations in blue and red, re¬ 
spectively. 

It is important to determine the ranges in effective ra¬ 
dius and magnitude over which the NGVS KINGPHOT 
measurements are reliable. We have addressed this ques¬ 
tion in several ways. Eirst, we carried out a large num¬ 
ber of simulations in which artificial sources were gen¬ 
erated using PSE-convolved King models plus a realistic 
amount of noise. The simulated objects were then added 
randomly to NGVS images and their sizes measured 
using KINGPHOT. The artificial objects had ^-band 
magnitudes in the range 20 to 23 mag, half-light radii 

These are equivalent to aperture radii of 0728, 0737, 0747, 
0756, 0765, 0775, 175 and 370. In physical units, these radii are 
22, 30, 37, 45, 52, 60, 120 and 240 pc at the distance of Virgo. 
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Fig. 1.— Aperture correction map (from 8-pixel diameter aper¬ 
tures to 16-pixel diameter apertures) made by stacking all individ¬ 
ual NGVS fields. The four panels show results for the u* (upper-left 
panel), g (upper-right panel), i (lower-left panel) and 2 ; (lower-right 
panel) bands. Color coding shows the deviations in magnitudes 
from the mean correction across each MegaCam pointing. 


between 1 and 100 pc, and concentration parameters, 
c = log;LoU/Fc, between 0.5 and 2.3^^. Figure 2 com¬ 
pares the actual and measured cluster sizes from these 
simulations; the symbols show mean values of Arh/vh 
while the errorbars show the standard deviation at each 
Vh- Results are shown, from top to bottom, in bins of 
decreasing magnitude (^ = 20 to 23 mag). The KING- 
PHOT measurements are found to be in good agreement 
with the actual sizes over a wide range in radius, with 
the exception of the most compact objects which have 
input sizes of rh ^ 2 pc. 

The simulations show that there is a tendency to over¬ 
estimate the sizes of the smallest objects and underes¬ 
timate the sizes of the largest objects. The most reli¬ 
able KINGPHOT measurements are obtained for objects 
brighter than g ^ 21.5 mag and having half-light radii 
of Vh r\j 30 pc. Of course, there are good reasons to ex¬ 
pect that the measurements for the larger {vh ^ 50 pc) 
objects would show increased scatter. First, at fixed 
magnitude, larger objects always have lower mean sur¬ 
face brightness, and the lower S/N ratio will naturally 
lead to an enhanced scatter. Second, in our analysis 
the KINGPHOT fitting radius for the NGVS is 7 pixels. 
This corresponds to 105 pc at the distance of Virgo. 
Jordan et al. (2005) have shown that the KINGPHOT 
measurements begin to show a bias when ^ Ffit/2, 
which translates to > 50 pc for the NGVS. As noted 
above, UCDs are often defined to have sizes between 10 
and 100 pc, although the exact limits are arbitrary and 
several recent studies have found that most UCDs have 
10 < r/i < 30 pc (Brodie et al. 2011; Chiboucas et al. 
2011; Strader et al. 2011; Penny et al. 2012). Our choice 


Here rt and rc are tidal and core radius, respectively. 
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Fig. 3.— A comparison of half-light radii measured in the g and 
i bands for objects brighter than g = 21.5. The dashed line shows 
the one-to-one relation, while the dotted lines show half-light radii 
of 5, 10, 50 and 100 pc. 

of fitting radius thus seems appropriate for the major¬ 
ity of known UCDs. We conclude from the simulations 
that the KINGPHOT measurements should be reliable 
for UCDs and GCs brighter than g 21.5 and having 
half-light radii in the range ^10-100 pc. 

A second check on the reliability of the KINGPHOT 
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Fig. 4. — Comparison of ^f-band half-light radii measured in 
the overlap regions of the four NGVS fields surrounding M87. All 
objects shown in this plot have g < 21.5. The dashed line shows 
the one-to-one relation while the dotted lines shows a half-light 
radius of = 10 pc. 

sizes is possible using the independent g- and i-band mea¬ 
surements. We compare the measurements in these two 
bands in Figure 3, now showing only those objects with 
g < 21.5 mag (see above). The dotted lines show sizes 
of 5, 10, 50 and 100 pc. The half-light radii are in good 
agreement over a wide range in radius — roughly from 
5 pc to 100 pc. A third check on the NGVS size mea¬ 
surements is possible using regions of overlap between ad¬ 
jacent NGVS fields, where some sources can have their 
sizes measured in completely independent stacks. Fig¬ 
ure 4 shows this comparison for the overlap regions in 
the four M87 fields. The agreement is once again very 
good, apart from the most compact objects. In fact, this 
comparison is likely to be overly pessimistic since the 
measurements in this case are made entirely at the edges 
of individual MegaCam fields where the image quality is 
poorest. 

Finally, we can compare our sizes with those measured 
from HST images. The blue squares in the upper panel of 
Figure 5 show the results of this comparison for GGs and 
UGDs from the AGSVGS (Jordan et al. 2005; Hasegan 
et al. 2005). Note that the NGVS sizes shown here refer 
to weighted mean sizes measured in the g and i bands. 
The agreement is quite good for objects with Vh ^ 10 pc, 
although it is worth bearing in mind that the AGSVGS 
sizes were also made with KINGPHOT. The cyan squares 
in this figure show objects from Strader et al. (2011), 
who collected HST images and measured half-light radii 
for spectroscopically confirmed GGs and UGDs using the 
ISHAPE program of Larsen (1999). For the smallest ob¬ 
jects {vh ^ 5 pc), the NGVS sizes are clearly overesti¬ 
mated, as expected from the simulations shown in Fig¬ 
ure 2. For larger objects, the NGVS measurements gen¬ 
erally compare favorably with the HST values, although 
there is a ^ 25% downward offset for Vh ^ 10 pc relative 
to the values of Strader et al. (2011). The explanation 
of this offset is unclear, and while it may be tempting to 
attribute the offset to the use of different algorithms (i.e.. 
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Fig. 5.— (Upper Panel) Comparison of CC and UCD half-light 
radii measured with the NCVS to those measured from HST-based 
imaging. The NCVS errorbars show the weighted mean KING- 
PHOT uncertainties in the g and i bands. The blue symbols show 
objects with published sizes from the ACSVCS (Jordan et al. 2005; 
Ha§egan et al. 2005; Jordan et al. 2009) while the cyan symbols 
show HST measurements from Strader et al. (2011) for CCs and 
UCDs in the M87 region. (Lower Panel) Comparison of NCVS 
half-light radii measured in the g and i bands for our final sample 
of UCD candidates in the M87, M49 and MOO regions. See §2.3 
and §2.4 for details on the selection procedure. 




Fig. 6 . — Comparison of the total UCD magnitude derived 
from KINGPHOT versus that from aperture-corrected, 16-pixel- 
diameter aperture photometry. The upper and lower panels show 
results for the g and i bands, respectively. The objects shown here 
are the final sample of UCDs identified in the M87 region (see §2.6 
for details). 


KINGPHOT vs. ISHAPE), previous studies have shown 
the codes to yield fully compatible results (e.g., Peng 
et al. 2009), at least for GG-sized sources. In any case, 
our aim in this study is not to measure UCD sizes at a 
level of precision that rivals HST, but simply to identify 
any UCDs with sizes in the range ^ 10 to 100 pc. 

After some experimentation, candidate UCDs identi¬ 
fied on the basis of their magnitudes and colors (§2.4) 













NGVS. X. Ultra-Compact Dwarfs in M87, M49 and M60 


7 


were considered to have sizes appropriate for bonafide 
UCDs if three conditions were satisfied: (1) the mea¬ 
surement errors on Vh were less than 1.5 pc, in both the 
g and i bands; (2) the fractional difference in the size 
measurements was \rh,g — Th,i\/{rh) < 0.5; and (3) the 
weighted mean half-light radius was (vh) > 11 pc. Addi¬ 
tionally, we limit our selection to sources brighter than 
g = 21.5 mag. Taken together, these criteria ensure that 
objects with uncertain sizes, spurious sources, and com¬ 
pact objects that may suffer from biased size measure¬ 
ments, are eliminated from our analysis. The lower panel 
of Figure 5 compares the g and i-band sizes for candidate 
UCDs identified in our three program regions. 

Before proceeding, we pause to explain the photomet¬ 
ric measurements used in this analysis. For each UCD 
candidate, we have u*giz aperture magnitudes (and, in 
the M87 region, grizKg)^ measured in a series of ever 
larger apertures as described above. We also have PSF- 
convolved King-model fits from KINCPHOT in the g and 
i bands. In principle, these latter values should provide 
the best photometry for the UCDs. However, there are 
two limitations with their use: (1) the KINGPHOT mea¬ 
surements are available only in two of the six bands; and 
(2) analogous measurements exist only for the bright¬ 
est {g < 21.5 mag) objects. For these reasons, in this 
paper we generally rely on the corrected, 16-pixel diam¬ 
eter aperture magnitudes when comparing the UCD and 
GC properties. In addition, aperture colors are most 
reliable when using relatively small apertures, so colors 
are based on measurements within an 8-pixel diameter 
aperture. Figure 6 compares the 16-pixel aperture and 
KINGPHOT magnitudes for our final sample of UCDs 
(see §2.6) in the M87 region, demonstrating that there is 
good agreement between the two sets of measurements. 


2.4. Sample Selection 

To select a sample of UCDs, we must rely on a com¬ 
bination of parameters: magnitudes, half-light radii, col¬ 
ors (specifically, their location within the color-color di¬ 
agram or diagrams), and mean effective surface bright¬ 
ness. The magnitude range adopted for the selection of 
both UCDs and GCs in our study is 18.5 < g < 21.5 
mag (—12.7 < Mg < —9.7). The bright limit is set by 
saturation in the NGVS long exposures while the faint 
limit reflects the minimum S/N needed to measure ac¬ 
curate sizes using KINGPHOT (see §2.3). Note that we 
select both UCD and GC candidates from their location 
in color-color diagrams (§2.4.1 and 2.4.2) and identify 
UCDs on the basis of their larger sizes (11 < < 100 pc, 

§2.3) and lower mean surface brightness at fixed luminos¬ 
ity (§2.4.2). 

In this work, we aim to study the properties 
of UCDs around three luminous early-type galaxies 
in the Virgo cluster: M87 (=NGC4486, VCC1316), 
M49 (=NGC4472, VCC1226) and M60 (=NGC4649, 
VCC1978). Some basic properties of these three galaxies 
are listed in Table 1, including coordinates, magnitudes, 
colors, effective radii, mean effective surface brightness, 
concentration cs 2 (defined as the ratio of the radii con¬ 
taining 80% and 20% of the luminosity derived from 
a curve-of-growth analysis), stellar luminosity, stellar 
mass. X-ray gas mass, dark matter mass, VCC morphol- 



(M87, M49 and M60) whose UCD populations are examined in 
this paper. Red squares show individual NGVS pointings. VCC 
galaxies that are confirmed or likely members of the cluster are 
shown as black circles, with symbol size proportional to luminosity. 
The solid and dotted lines show the respective boundaries of the 
NGVS and VCC surveys. The large circles denote the virial radii 
of the Virgo A and B subclusters (see Ferrarese et al. 2012). 

ogy, NGVS classification type, and other names. 

Although these are the three brightest galaxies in the 
cluster^^, they occupy very different environments, as 
shown in Figure 7. As is well known, Virgo is not a dy¬ 
namically relaxed system but consists of several distinct 
subclusters (e.g., Binggeli et al. 1987, 1993; Schindler 
et al. 1999; Mei et al. 2007). Subcluster A, which is 
centered on M87, is the dominant component of Virgo. 
Subcluster B, which is centered on M49 and located ^4° 
to the south of M87, is considerably smaller and more 
compact than subcluster A (Binggeli et al. 1993). The 
third, less conspicuous, subcluster (C) is centered on a 
small group of galaxies containing M59 and M60 located 
^3° to the east of M87 (Binggeli et al. 1987).^^ Simi¬ 
larly, X-ray imaging of the cluster shows some dramatic 
differences in the properties of the hot gas in the vicin¬ 
ity of these galaxies (Bohringer et al. 1994). Located at 
the center of subcluster A, M87 is unsurprisingly embed¬ 
ded in a vast, massive reservoir of X-ray-emitting gas. 
M49 contains much less (^ 1/40) hot gas than M87, but 
substantially more than M60 (4-5x). 

These three galaxies also reside in quite different envi¬ 
ronments in terms of their surrounding diffuse light. Lo¬ 
cated at the cluster center, M87 possesses an extended 
diffuse envelope and is embedded in a web of intracluster 
light (Mihos et al. 2009; Castro-Rodriguez et al. 2009; 
Rudick et al. 2010). However, this web is highly concen¬ 
trated in the central regions of Virgo (Castro-Rodriguez 
et al. 2009; Boselli et al. 2014, Mihos et al. 2015, in 
preparation). Rather than being embedded in an ex- 

The blue luminosities of the three galaxies are quite similar, 
differing from their mean luminosity by only < 25%. 

A fourth subcluster is probably associated with M86 
(VCC763), located ~ 1?5 WNW of M87 (e.g., Schindler et al. 1999; 
Binggeli 1999; Mei et al. 2007). 
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TABLE 1 

Basic Data for M87, M49 and M60. 



Parameter 

(1) 

Units 

(2) 

M87 

(3) 

M49 

(4) 

M60 

(5) 

Reference'^ 

(6) 

a(J2000) 

deg 

187.7059365 

187.4448543 

190.9165323 

1 

(5(J2000) 

deg 

+12.3911224 

+8.0004889 

+11.5527000 

1 

9 


mag 

9.03 

8.71 

9.11 

1 

{(9- 

i)) 

mag 

1.15 

1.12 

1.10 

1 

Re 


arcsec 

96.3 

105.7 

76.0 

1 

Me 


mag arcsec”^ 

20.95 

19.67 

19.22 

1 

C82 



10.06 

9.85 

9.71 

1 

logio 



10.92 

11.04 

10.89 

1 

logio 

[M^/Mq] 


11.55 

11.63 

11.46 

1 

logio 

Mgas/^Q] 


13.28 

11.64 

-11 

3,4 

logic 

Mdm/^o] 


14.15 

13.94 

13.54 

3,4 

VGG Morphology 


E2/S0i(2) 

EO 

S0i(2) 

2 

NGVS Glassification Type 


E-Q 

E(sh st)-Q 

E(sh:)-Q 

1 

NGG 



4486 

4472 

4649 


VGG 



1316 

1226 

1978 

2 


't'Key to references: (1) NGVS ; (2) Binggeli et al. (1985); (3) Schindler et al. (1999); (4) Humphrey 
et al. (2006). All magnitudes and surface brightnesses have been extinction corrected and transposed 
to the SDSS photometric system. Luminosities and masses are calculated assuming the distance moduli 
listed in Mei et al. (2007). 


tended envelope of diffuse light, M49 shows tidal features 
scattered througout its halo (Janowiecki et al. 2010; Ar- 
rigoni Battaia et al. 2012; Ferrarese et al. 2012), while 
MOO shows diffuse light associated with M59 (Yan et al. 
2008). These differences in the prominence and morphol¬ 
ogy the diffuse light surrounding each galaxy are likely 
related to the different dynamical histories of the three 
regions. 

2.4.1. u^iKs Selection 

In the NGVS pilot program region (the four helds in 
the immediate vicinity of M87: NGVS+0+0, NGVS+0+1, 
NGVS-1+0 and NGVS-l+l), we have imaging in ah hve op¬ 
tical bands, u^griz (Ferrarese et al. 2012), plus LCg-band 
imaging from the NGVS-IR (Munoz et al. 2014). In this 
region, we are thus able to select UCD and GC candi¬ 
dates directly from the (i^* — i)-{i — Kg) color-color dia¬ 
gram (= u*iKs) which was shown by Munoz et al. (2014) 
to be remarkably effective in separating UCDs and GGs 
from foreground stars and background galaxies. Indeed, 
Figure 8 shows that the sources in the M87 region — 
many of which have been observed spectroscopically (see 
Zhu et al. 2014; Zhang et al. 2015 — can be neatly sepa¬ 
rated into three main groups in the u*iKs diagram: GGs 
and UGDs, background galaxies, and foreground stars. 
The polygon in this hgure shows our adopted selection 
region for UGDs and GGs, a region that includes ah of 
the spectroscopically conhrmed UGDs and most of the 
conhrmed GGs in the vicinity of M87 (Hanes et al. 2001; 
Strader et al. 2011; Zhang et al. 2015). 

2.4.2. u'^gz Selection 

Unfortunately, of the 117 MegaCam helds that make 
up the NGVS, a LCg-band photometric catalog is avail¬ 
able at this time in only the four pilot program helds. For 
the remaining helds, including the M49 and M60 regions, 
we must rely on the (rt* — g)-{g — z) color-color diagram 
(= u^gz) to identify UGD and GG candidates. While 
there is certainly a separation of UGDs and GGs from 
foreground stars and background galaxies in this color- 



Fig. 8.— Distribution of sources in the M87 region in the (n* -i)- 
{i — Ks) diagram (gray symbols). Blue squares, cyan crosses, green 
stars and orange triangles show the location of spectroscopically 
confirmed UCDs, GCs, stars, and background galaxies, respec¬ 
tively. 

color space, it is not as dramatic as in the case of u*iKs 
(Munoz et al. 2014). As Figure 9 demonstrates, the sam¬ 
ple of selected UGD and GG candidates (blue and 
cyan dots) shows some overlap with the regions occupied 
by spectroscopically conhrmed stars and galaxies (green 
and orange dots). As the previous hgure, the polygon 
in Figure 9 shows our adopted selection region for UGDs 
and GCs. 

To better understand contamination arising from the 
selection of UCDs and GCs in the u'^gz diagram, we 
show in Figure 10 the distribution of 4t*^^-selected can¬ 
didates from the M87 region in the u^iKg diagram (red 
dots). For comparison, the complete sample of sources 
in this region is shown by the faint gray dots. Clearly, 
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Fig. 9. — Distribution of sources in the M87 region in the (n* -g)- 
(g — z) diagram (gray symbols). Blue squares, cyan crosses, green 
stars and orange triangles show the location of spectroscopically 
confirmed UCDs, GCs, stars, and background galaxies, respec¬ 
tively. 



Fig. 10. — Distribution of gz-selected sources from Figure 9 in 
the u*iKs diagram. Red and blue symbols denote the combined 
sample of UCD and GC candidates before and after a second se¬ 
lection on mean effective surface brightness. 

the level of contamination is non-negligible, with many 
of the ix*^^-selected UCD and CCs candidates located in 
the regions occupied by stars or galaxies. 

To reduce the contamination in this u*gz-selected sam¬ 
ple, we impose an additional selection on mean effective 
surface brightness, (/i^)e- In the upper panel of Fig¬ 
ure 11, the spectroscopically confirmed UCDs and CCs 
are found to occupy a relatively narrow range in sur¬ 
face brightness (dashed histogram). The lower panel of 
Figure 11 shows mean effective surface brightness plot¬ 
ted against {g — i) color for all u*gz- and selected 



Fig. 11. — (Upper Panel). Histogram of mean effective surface 
brightness, {gg)e, for sources in the M87 region. The solid his¬ 
togram shows the distribution for color-selected (both u^iKg- and 

gz-selected) UGDs and GG candidates, while the dashed his¬ 
togram shows the distribution for the subset of spectroscopically 
confirmed UGD and GGs. The dotted lines show our adopted sur¬ 
face brightness limits for the selection of UGDs. (Lower Panel). 
Dependence of (g — i) color on mean effective surface brightness. 
Gray dots show all u*gz-selected UGD and GG candidates. Spec¬ 
troscopically confirmed UGDs, GGs and galaxies are shown by blue 
squares, cyan crosses and orange triangles, respectively. 

sources from the M87 region (gray and black dots, re¬ 
spectively). Spectroscopically confirmed UCDs, CCs and 
galaxies are shown by the blue, cyan and orange sym¬ 
bols, as indicated in the legend. The irregular pentagon 
(denoted by the dotted lines) shows the region of (/i^)e- 
{g—i) space used to identify UCDs and CCs; the diagonal 
boundary in the upper right corner was used to minimize 
contamination from background galaxies, which tend to 
have red colors and low surface brightness. While a selec¬ 
tion on surface brightness does mean that a small number 
of (high surface brightness) CCs will be missed, all but 
four of the spectroscopically confirmed UCDs are identi¬ 
fied using this approach. 

Figure 10 plots the u*gz and surface brightness- 
selected objects in the u*iKs diagram as blue dots. Al¬ 
though a handful of objects are still found in the regions 
of the diagram occupied by stars or galaxies, the sample 
is clearly much cleaner than before. We notice that the 
u^iKs-seiection does introduce many background galax¬ 
ies, as seen in the lower panel of Figure 11 (the cloud 
of black points in the upper right corner). Therefore, 
we proceed by selecting UCD candidates relying on a 
combination of u*iKs and {iJig)e-{g — i) diagram for the 
M87 region, and a combination of u*gz and {fig)e-{g — i) 
diagrams for M49 and M60 regions. 

2.4.3. Visual Inspection of UCD Candidates 

After the color-color selection, the UCD candidates in 
all three regions (M87, M49 and M60) were separated 
from CCs on the basis of their half-light radii. Oper¬ 
ationally, the UCDs were defined as those objects with 
r/i > 11 pc (§2.3; see also Drinkwater et al. 2004; Will- 
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Fig. 12.— UCD candidates 1-25 in the M87 region selected on 
the basis of magnitude, color and half-light radius (see Tables 2 
and 3). Most of these candidates are class 1 objects, although 
some objects are known nucleated dwarf galaxies (7, 9, 12) or 
likely background galaxies (15). Each panel measures 60x60 pixels 
(1172x1172). North is up and east to the left. 

man & Strader 2012; Briins & Kroupa 2012; Penny et al. 
2012) and Mg < —9.7 mag. Up to this point, UCD se¬ 
lection was performed in a completely objective manner. 
However, a final, somewhat subjective step — a visual in¬ 
spection of each candidate — is required to eliminate ob¬ 
vious contaminants or artifacts. Such contaminants were 
usually found to be nucleated dwarf galaxies or compact 
background galaxies. After a visual inspection, each can¬ 
didate was assigned a classification code to identify it as 
a UCD or an another type of object: (1) = probable 
UCD; (2) = dwarf nucleus; (3) = background galaxy; 
(4) = other type of contaminant (usually a candidate 
UCD with uncertain properties due to its being located 
close to a bright star, stellar halo or CCD bleed trail). 

Figures 12-16 present ^-band mosaics for our complete, 
objectively selected sample of candidate UCDs in the 
M87 region. Each panel has been labelled with both 
our class parameter and the object identification num¬ 
ber from Tables 2 and 3. These 127 objects make up our 
sample of UCD candidates that satisfy the joint selection 
criteria on magnitude, color and size described above. 

2.4.4. Constraints on Membership and Contamination 
from Speetroseopy 

While the selection criteria described above can pro¬ 
vide a very clean sample of likely UCDs from NGVS 
imaging alone, there is no substitute for spectroscopic 
redshifts in distinguishing UCDs from foreground stars 
or background galaxies. 

Of course, there is no single spectroscopic survey of the 
Virgo cluster that is complete to the depth of the UCD 
samples considered here, but there nevertheless exists a 
huge amount of spectroscopic data for bright, compact 
sources in the direction of Virgo, especially in the clus¬ 
ter core. Such studies include early surveys carried out 
with MMT (Huchra & Brodie 1987) and Palomar (Mould 


Fig. 13. — Same as Figure 12 except for objects 26-50. 
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Fig. 14. — Same as Figure 12 except for objects 51-75. 


et al. 1990), as well as subsequent programs with Keck 
(Cohen & Ryzhov 1997; Cohen 2000; Strader et al. 2011), 
CFHT (Hanes et al. 2001) and AAT (Jones et al. 2006). 
Beyond the immediate vicinity of M87, the Sloan Digital 
Sky Survey (SDSS; York et al. 2000) observed many tar¬ 
gets over the full extent of the cluster, although these 
were mainly Milky Way stars, QSOs and background 
galaxies. Most recently, a number of coordinated NGVS 
spectroscopic programs have been undertaken with the 
MMT, AAT and other facilities (see, e.g., Zhang et al. 
2015). Thus, there exists a rather extensive spectroscopic 
database that can be compared to our UGD catalog. 

Among the sample of 127 UGD candidates found in 
the M87 region (before visual inspection), 85 have mea¬ 
sured velocities which are listed in Table 3. Using these 
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Fig. 15. — Same as Figure 12 except for objects 76-100. 
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Fig. 16. — Same as Figure 12 except for objects 101-127. 


85 objects, we will estimate the contaminants in both 
the u^iKg- and ix*^z-selected UCD samples. Of these, 70 
and 73 candidates, respectively, were deemed to be likely 
UCDs on the basis of u'^iKs and u'^gz diagrams, and were 
classified as class = 1 systems. Four of 70 objects in 
the IX*-selected sample and eight of 73 objects in the 
XX*selected sample have velocities {vr > 3500 km/s) 
that show them to be background galaxies. We would 
expect ^ 6% of xx*xi^s-selected candidates and ^ 11% of 
XX* selected candidates remaining after visual inspec¬ 
tion to be misclassified background galaxies. 


As noted in §1, a companion paper in the NGVS series 
presents a kinematic study of UCDs in the M87 region 
(Zhang et al. 2015). Since their sample is a little dif¬ 
ferent from ours, we pause to compare the slightly dif¬ 
ferent selection strategies used in the two studies. The 
Zhang et al. (2015) analysis is based on radial velocities 
for 97 UCDs which were identified using a combination 
of NGVS and HST imaging. Among this sample, 86 
UCDs with radial velocity measurements were selected 
from the NGVS to have > 11 pc and g < 21.5. Eleven 
more UCDs^^ with measured radial velocities were also 
included in the Zhang et al. 2015 analysis. Based on 
earlier, accurate HST size measurements, some of these 
objects have sizes in the range 9.5 < r/,, < 11 pc, which is 
just below the adopted NGVS cutoff in r/j,, while the oth¬ 
ers are larger than r/j, = 11 pc but fainter than g = 21.5 
mag. In short, Zhang et al. (2015) adopted slightly more 
flexible selection criteria in order to explore the kinemat¬ 
ics of the M87 UCD system using the largest possible 
radial velocity sample. Our sample selection relies en¬ 
tirely on NGVS imaging in an effort to maximize spatial 
completeness and sample homogeneity. 

2.5. Control Fields and Estimation of Contamination 

The NGVS includes four background (BG) fields lo¬ 
cated far from the cluster center. Each field is offset 
from M87 by ^16°, or roughly three virial radii.Be¬ 
cause these four fields were observed in the same filters as 
the other survey fields (rCgriz)^ and under typical NGVS 
observing conditions with identical exposure times and 
dithering strategies, we can use these fields to estimate 
the level of contamination in our UCD samples selected 
from xx*^z imaging. 

In all, the background fields contain a total of 10 322 
objects between our UCD magnitude selection limits of 
18.5 < g < 21.5. Of these, 1145 objects are addi¬ 
tionally located in the region of the xx*^z diagram oc¬ 
cupied by UCDs and GCs. The selection on surface 
brightness described in §2.4.2 is quite effective at re¬ 
moving galaxies, however, and leaves a total of only 174 
sources. Among these 174 sources, 19 were found to 
satisfy the three size conditions described in §2.3. Af¬ 
ter a visual inspection, only 9 “UCDs” were left in the 
sample, a slightly higher number than the four spuri¬ 
ous sources identified based on spectroscopy (§2.4.4) of 
xx*iArs-selected sources. While the field-to-field variance 
may be significant (there were 2, 1, 4 and 2 contaminants 
found in BGl, BG2, BG3 and BG4 fields, respectively), 
the number density of contaminants in the xx*^z-selected 
sample is E(UCD)u*gz = 2.2 ± 0.7 deg“^. 

Unfortunately, there is no ATg-band imaging for these 
four background fields. To estimate the contamination in 
our XX* xATg-selected sample, we use the CEHTLS-D2 field 
as a proxy background field for the NGVS. CEHTLS-D2 
field is one of the four CEHTLS Deep Survey pointings 
observed in u*griz filters. In this region, there is also 
a deep Kg imaging from UltraVISTA (McCracken et al. 
2012) which is an ultra-deep, near-infrared survey us¬ 
ing VISTA telescope. We find 2731 objects in this field 

28 of these 11 objects were not detected in the NGVS due 

to their proximity to bright stars. 

We note that none of these fields contains a known galaxy 
cluster with z < 0.1, although a galaxy group, located at 2 : ~ 0.02 
(Crook et al. 2007), does fall within the BG4 field. 
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in the magnitude range 18.5 < g < 21.5 mag; reassur¬ 
ingly, this is quite similar to the mean number of 2580 
objects per NGVS background field in this range. A 
total of 36 CFHTLS objects fall within the UCD/GC re¬ 
gion of the u^iKs diagram, only one of which remains 
after selecting on size and surface brightness. Thus, 
the i4*iArs-selected UGD sample is found to be excep¬ 
tionally clean, with a contaminant surface density of 
S(UGD)u*iK3 ^ 0.9 ± 0.9 deg-^ 

2.6. Catalog of UCD Candidates 

The complete sample of 127 UGD candidates in the 
M87 region is presented in Tables 2 and 3. Recall that 
this sample was selected entirely on the basis of mag¬ 
nitude, color (specifically, location in the and/or 

u'^gz diagrams), half-light radius and surface brightness, 
as described in the preceding sections. For complete¬ 
ness, we list all candidates that were selected before vi¬ 
sual inspection and including even those objects which 
are known from spectroscopy to be background galaxies. 
A total of 16 of these 127 UGD candidates were subse¬ 
quently rejected following a careful examination of the 
NGVS images (i.e., the objects with class parameters 
of 2, 3 or 4 in Table 3). 

Information on the M87 UGD candidates is divided 
into two tables. Table 2 gives, from left to right, each 
object’s identification number, NGVS identifier, right as¬ 
cension, declination, extinction (Schlegel et al. 1998), g- 
band magnitude measured in an aperture with a diam¬ 
eter of 16 pixels and corrected to the SDSS PSF mag¬ 
nitude system, and a variety of 8-pixel diameter aper¬ 
ture colors. Note that these are the observed magni¬ 
tudes, prior to any correction for extinction. The final 
two columns of Table 2 give the parameter Aenv, defined 
in §4, followed by our final UGD classification: “1” if 
the candidate is a probable UGD after considering all 
possible diagnostics (magnitude, color, size, morphology 
and/or spectroscopy) and “0” otherwise. A total of 92 
objects have UGD classifications of 1. 

Table 3 records the object identification number, 
NGVS identifier, g- and i-band KINGPHOT magnitudes 
and radii, followed by previous measurements from the 
AGSVGS (Jordan et al. 2005; Hasegan et al. 2005) or 
from Strader et al. (2011). The next four columns give 
the measured radial velocity, if available, from Strader 
et al. (2011), MMT or AAT (Zhang et al. 2015), and 
SDSS. The NGVS class parameter is given in the next 
column, followed by two flags to indicate whether the 
object falls inside (1) or outside (0) the GG/UGD se¬ 
lection regions in the u'iKs and u'gz diagrams, re¬ 
spectively, and the UGD classification parameter de¬ 
fined above. The final column gives alternative names 
from the literature (Binggeli et al. 1985; Hanes et al. 
2001; Jones et al. 2006; Firth et al. 2008; Harris 2009a; 
Zhang et al. 2015). One other object (ID = 7 = NGVS- 
J122643.32+121743.9) was identified as a faint, nucleated 
galaxy in our search for low-mass galaxies in the Virgo 
cluster core. 

In the M87 region, 106 of the 127 candidates listed 
in Tables 2 and 3 were selected based on magnitude, 
location in R'iKs diagram, size and surface brightness 
{uiKs = 1 in Table 3). Among the 106 Kg-selected 
candidates, 96 were classified as class = 1 objects. Ra¬ 
dial velocities reveal four of these objects to be back¬ 


ground galaxies (see §2.4.4), bringing our cleaned sample 
of UGDs in the M87 region to 92. Our cleaned sample 
of 92 UGD candidates is much larger than any previous 
studied UGD samples in this region (e.g. 34 UGDs in 
Brodie et al. 2011 and Strader et al. 2011). In addition, 
the homogeneity is another advantage of our sample. In 
addition, we also have a homogeneous sample GG and 
dE,N (see §2.7 and §2.8), which allows us to compare the 
properties between UGDs, GGs and dwarf nuclei. 

As described in §2.4.2, we also select UGD candi¬ 
dates based on the R'gz diagram in the M87, M49 and 
MOO regions. Based on magnitude, location in R'gz 
diagram, size and surface brightness, we find 125, 50 
and 29 UGD candidates, and 110, 40 and 23 are clas¬ 
sified as class = 1 in these three regions. 8 and 12 
candidates in the M87 and M49 regions have velocity 
Vr > 3500 km/s, which we deem to be compact, back¬ 
ground galaxies. For the MOO region, only two candi¬ 
dates have velocity measurements and both of them are 
Virgo members {vr < 3500 km/s). Finally, the cleaned 
i4*^z-selected samples have 102, 28 and 23 confirmed or 
candidate UGDs in these three regions. Table 4 and 5 
show the photometric and structure properties of M49 
UGD candidates, while Table 6 and 7 show the proper¬ 
ties of M60 UGD candidates. 

The properties of our final UGD samples are summa¬ 
rized in Table 8, including the region geometry, coordi¬ 
nates, total area, and the number and surface density of 
the objects within each region. Note that the numbers 
given in this table refer to the total numbers of candi¬ 
dates identified using the different color-color diagrams, 
with no selection on redshift. 

It is difficult to find all the UGDs within specific mag¬ 
nitude, color, size and surface brightness range. Based 
on our selection criteria, we cannot find bright UGDs 
{g < 18.5) that are saturated in NGVS long exposure 
images. We may also miss UGDs that are located close 
to saturated stars or galaxy centers. In addition, the 
half-light radii of UGDs might increase with galactocen- 
tric distance (Pfeffer & Baumgardt 2013), which can 
introduce a location-dependent incompleteness via our 
adapted constant size cut Vh ^ 11 pc. Finally, uncer¬ 
tainties on the size and surface brightness measurements 
could influence the sample completeness. We know, 
for instance, that we miss four UGDs due to our sur¬ 
face brightness cut (see Figure 11). A fixed concentra¬ 
tion c = 1.5 is assumed (the same as in Jordan et al. 
2005) when we measure the half-light radius using KING- 
PHOT . Thus, the radii of UGDs with larger concen¬ 
trated (c > 1.5) are underestimated, and vice versa. In 
this case, more compact UGDs (c > 1.5) may be missed 
via our size criteria. The PSF variation across the survey 
area could also influence our size measurements. How¬ 
ever, the variation is small (see Figure 8 of Ferrarese 
et al. 2012) and we measured the half-light radii by fit¬ 
ting PSF-convolved King (1996) models. So the effect of 
PSF variations on our UGD sample selection is small. 

2.7. Globular Cluster Sample 

As seen in Figure 8, the GGs and UGDs are clearly 
separated from foreground stars and background galax¬ 
ies. In the M87 region, it is reasonable easy to select 
GG candidates brighter than ^ = 24 mag in the u'^iRg 
diagram (Munoz et al. 2014). We did not adopt other 
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criteria like surface brightness and size since it is diffi¬ 
cult to measure the sizes of objects down to ^ = 24 mag, 
and most GCs are unresolved. We eliminate the UCD 
candidates described above and any other known back¬ 
ground galaxies and foreground stars. Finally, we select 
9497 GC candidates which are brighter than ^ = 24 mag. 

Around M49 and M60, where we do not have Kg 
photometry, the GCs are selected based on u^gz color- 
color diagram. The u^gz diagram is not as optimal as 
u^iKg to separate GCs from stars and background galax¬ 
ies, especially for the faint objects. Therefore, we only 
select bright GCs in the M49 and M60 regions. The 
adopted magnitude criteria is the same as UCD selec¬ 
tion, 18.5 < ^ < 21.5 mag. Finally, 307 and 273 bright 
GCs are found in the M49 and M60 regions respectively. 
In this paper, we only use the M49 and M60 bright GCs 
to compare with UCDs, and do not study the proper¬ 
ties of GCs that need a complete sample, e.g. luminosity 
function. In this case, the magnitude completeness of 
our GC sample is not very important. 

2.8. Nucleated Galaxies in the M87 and M49 Region 

A leading theory for the origin of UCDs is tidal strip¬ 
ping, or “threshing”, of nucleated dwarf galaxies. In this 
scenario, the outer envelope of a low-mass satellite is re¬ 
moved by the strong tidal field of a cluster or galaxy, leav¬ 
ing behind a nucleus that survives because of its compact 
nature (e.g., Bekki et al. 2001, 2003; Drinkwater et al. 
2003; Pfeffer & Baumgardt 2013; see also Bassino et al. 
1994 for an early exploration of this same basic idea). 
It is therefore of interest to compare the properties of 
UCDs and nuclei in the NGVS. Accordingly, we pause 
to describe the construction of a sample of nuclei for the 
M87 region that will be used for such comparisons. 

The M87 region contains a total of 89 certain or possi¬ 
ble nucleated galaxies from either Binggeli et al. (1985) or 
the NGVS (Ferrarese et ah, in preparation). Nucleated 
galaxies are identified based on visual inspection of the 
images as well as parametric fits to the surface brightness 
profiles. Galaxies are deemed to be “certainly nucleated” 
when the images show clear evidence of a central lumi¬ 
nosity excess on the scale of the PSF, and the measured 
surface brightness profiles lie above over the inward ex¬ 
trapolation of the Sersic law best fitting the profile be¬ 
yond a few arcseconds. Galaxies for which either the 
visual inspection or the profile fits are inconclusive are 
classified as “possibly nucleated”; all others are classified 
as non nucleated (see also Gote et al. 2006 and Turner 
et al. 2012). Among the 89 certain or possible nucleated 
galaxies, nine objects were identified as UGD candidates 
in our automated selection process. While this provides 
direct evidence that at least some nuclei are robustly 
similar to the UGDs studied here, it is important to un¬ 
derstand why the remaining 80 objects do not appear in 
our UGD sample. First, nine of the 89 galaxies were not 
detected by SExtractor (due either to low surface bright¬ 
ness or projection close to a bright source), while 16 are 
brighter than g = 18.5 mag and 23 more are fainter than 
g = 21.5 mag, meaning that 41 nuclei remain after our 
selection on magnitude. Among these 41 nuclei, 6 of 
them fall outside the selection region in the u'^iKg dia¬ 
gram, 18 were rejected according to their surface bright¬ 
ness (most of them have lower surface brightness than 
UGDs) and another 8 were excluded when applying our 


rather stringent selection on size, leaving an objectively- 
selected sample of 9 nuclei finally. 

In what follows, we shall use both the objectively- 
selected sample of nine galaxy nuclei and the full sample 
of 64 nuclei^^ when comparing the UGDs and nuclei in 
the M87 region. In the M49 region, we find 22 galactic 
nuclei that have magnitude 18.5 < g < 21.5 mag. Two 
of them recovered as UGDs by our selection criteria. 

3. RESULTS 

The goals of our study are an improved understanding 
of the origin of UGDs and a comparison of their proper¬ 
ties to those of other stellar systems in the cluster, such 
as galaxies and GGs. We therefore begin by examining 
the bulk properties of the UGDs in the M87, M49 and 
M60 regions in order to understand the role environment 
may play in their formation and evolution. In a similar 
vein, we compare the properties of UGDs in these fields 
to those of GGs to determine the extent to which GGs 
and UGDs may share common formation and evolution¬ 
ary histories. Since the sample of UGDs in the M60 re¬ 
gion is smaller, some results in this section focus on just 
the M87 and M49 regions, e.g., color distribution, color 
gradient, color magnitude relation, and the connection 
to the Intracluster light. 

3.1. Population Richness 

Because the UGD samples identified in the three re¬ 
gions were selected using different approaches (i.e., the 
primary diagnostic is position in the u^'iKs diagram for 
M87, versus location in u'^gz for M49 and M60), we may 
expect some differences in the level of background con¬ 
tamination (see Table 3). The number of contaminants 
can be estimated robustly using either approach (§2.5), 
so we can straightforwardly compare the most basic pa¬ 
rameter for the UGD systems in these regions: the to¬ 
tal number of UGDs, Vucd- To ensure the most ho¬ 
mogenous estimates for Vucd, we focus on the results 
obtained when the i4*^2:-selection is used in all three re¬ 
gions. 

As described in §2.6, we have a total of 102, 28 and 
23 cleaned selected UGD candidates in the M87, 
M49 and M60 regions, respectively. To compare Vucd 
in the different regions, we count the number of UGD 
candidates within 2bRe of each galaxy^^. Finally, we 
find VucD ^ 78 ± 9, 22 ± 6 and 10 ± 4 within 25i?e of 
M87, M49 and MOO. Thus, it is immediately apparent 
that there are large region-to-region differences in pop¬ 
ulation size, despite the fact that the central galaxies 
have similar luminosities and stellar masses (Table 1). 
Indeed, in the MOO region, the detection of a UGD sys¬ 
tem is only significant at the ^2.5cr level. However, in 
the case of M49, and especially M87, the detection of a 
UGD system is unambiguous. 

Figure 17 shows the variation in A5 jcd with several 
fundamental parameters for the three regions of the 
cluster. Panels (a) through (f) show the dependence 
on galaxy ^-band luminosity, stellar mass, M*, X-ray 

These 64 objects refer to the sample of 80 nuclei that appear 
in the SExtractor catalog after discarding the 16 bright {g < 18.5 
mag) nuclei that are saturated in the deep NGVS images. 

We choose 25Re to make sure that the N\jcd count areas are 
in our survey regions. 
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Fig. 17.— Scaling relations for the number of UCD candidates, 
-^UCD, within 25Re of each galaxies. Panels (a) to (f) show the 
variation in A^ucd as a function of galaxy ^f-band luminosity, L*, 
stellar mass, M*, X-ray gas mass, Mgas, dark matter mass, Mdm? 
number of galaxies brighter than B = 18 within a radius of 25 Re, 
and the total number of globular clusters. The dashed lines show 
the best fit linear relations. 


gas mass, Mgas (Schindler et al. 1999; Humphrey et ah 
2006), dark matter mass, Mdm (Schindler et al. 1999; 
Humphrey et al. 2006), number of Virgo member galax¬ 
ies located within a radius of 2bRe, ^gab and the total 
number of GCs belonging to each galaxy, Nqq (Peng 
et al. 2008). As already noted, there is no correlation 
between the luminosity or stellar mass of the galaxy and 
the number of surrounding UCDs: i.e., M87 and MOO 
have nearly identical luminosities and stellar masses, but 
M87 contains a rich UCD system of ^80 objects, while 
MOO shows just a 2 . 5(7 enhancement in the number of 
UCD-like objects expected in a random field. Thus, it is 
clear that the MOO region is markedly deficient in UCDs 
compared to M87, or even M49, which again has a com¬ 
parable luminosity and stellar mass. 

On the other hand, the remaining panels in this figure 
show surprisingly tight correlations between VucD and 
the parameters describing mass of the host sub-cluster. 
The dashed lines show the best-fit linear relation in each 
case: 


log VuCD 


(0.380 ± 0.031) log Mgas - 3.14 ± 0.11 
(1.382 ± 0.321) log Mdm - 17.77 ± 1.20 
(1.581 ± 0.395) log Vgai - 0.70 ± 0.46 
(1.829 ± 0.092) log Vgc - 5.74 ± 0.18 


( 2 ) 

Evidently, the large differences in UCD population size 
and formation efficiency disappear when Vucd is com¬ 
pared to tracers of sub-cluster mass and/or richness. 
Normalized to the parameters such as dark matter mass 
and the number of neighboring galaxies, all three regions 
seem to have formed UCDs with comparable efficiency. 


• candidate globular cluster (g ^ 24.0) Ofull sample of dwarf nuclei 
■ spectroscopically confirmed UCD •dwarf nucleus in our sample 

□ candidate UCD (u'iK) + Virgo cluster member (B ^18) 



Fig. 18.— Spatial distribution of candidate UCDs (open squares) 
in the M87 region. Filled squares indicate UCDs that are spectro¬ 
scopically confirmed cluster members. CCs brighter than = 24 as 
shown as gray dots. Blue crosses indicate probable member galax¬ 
ies brighter than B = 18 from Binggeli et al. (1985). Red circles 
show nucleated dwarf galaxies in this region identified from the 
NGVS (see §2.8). Thin black lines show the ellipses that best fit 
the galaxy isophotes at [ig = 20, 21, ... 28 mag arcsec“^. The green 
curve shows the ellipse that best fit the global UCD distribution 
(see §3.2 for details). 


3.2. Spatial Distribution 

The spatial distribution of the confirmed and candi¬ 
date UCDs in each of these three regions are shown in 
Figures 18, 19 and 20. In each case, the distribution 
of UCDs (filled or open squares, depending on whether 
the object has been confirmed spectroscopically) is shown 
along with that of the GCs (gray dots) and Virgo member 
galaxies brighter than B = IS (blue crosses). For com¬ 
parison, the low-mass galaxies whose nuclei were identi¬ 
fied as possible UCDs in our automated selection pipeline 
and full sample of nucleated dwarf galaxies in this region 
(as described in §2.8) are also shown in Figure 18 with 
filled and open red circles, respectively. 

The ellipses in each figure show the best-fit galaxy 
isophotes at jig = 20,21, ... 28 mag arcsec“^. In the 
case of M87, where GCs can be selected with a high level 
of confidence from the u^'iKs diagram, we plot all GC 
candidates brighter than ^ = 24. For M49 and MOO, no 
LCg-band imaging is available so we plot only those GCs 
selected from the u^'gz diagram with g < 21.5 (see §2.7). 
For comparison. Figure 21 shows the spatial distribution 
of GC- and UCD-like objects in the four background 
fields (also selected from the u'^gz diagram). This last 
figure gives an indication of the level of contamination 
we might expect in the M49 and MOO regions, as well as 
the field-to-field variance. BG3 contains four UCD-like 
objects, whereas BG2 contains just one. 

The UCD candidates in Figure 20 are only weakly con¬ 
centrated to the center of MOO, and at least a few of the 
candidates may be better associated with M59, a mas- 
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candidate globular cluster (g<21 .5) o full sample of dwarf nuclei 
■ spectroscopically confirmed UCD • dwarf nucleus in our sample 
□ candidate UCD (u*gz) + Virgo cluster member (B<18) 



Fig. 19.— Same as Figure 18 except for the M49 region. 

sive early-type galaxy located about 25' 120 kpc) to 

the northwest. All in all, the number of UCD candidates 
in the M60 region appears only slightly larger than what 
one would expect in a random field (as noted in §3.1, we 
estimate a total UCD population of ^ 10±4 above back¬ 
ground) , although the existence of at least some bonafide 
UCDs in this galaxy is clear: e.g., the high-mass object, 
M60-UCD1, is projected just ^ 1.'5 from the center of 
the galaxy (Strader et al. 2013). 

The situation is quite different in M49 and, especially, 
M87, where there is a clear excess of confirmed or candi¬ 
date UCDs at the position of the central galaxies. In the 
latter case, the UCD sample is large enough to permit a 
more quantitative examination of its flattening and ori¬ 
entation. The dashed lines in Figure 22 show the radial 
variation in ellipticity, e, and position angle, 6>, for the 
isophotes of M87 according to Janowiecki et al. (2010). 
Although our sample of UCDs is too small to measure 
reliable profiles for these parameters, the mean values 
obtained using the UCDs within a radius of 2bRe 


(e)ucD — 0.34 
(^)ucD — 129° 


( 3 ) 


are in good agreement with those of the galaxy/ICL 
(intracluster light). Note that the UCD measurements 
shown in Figure 22 (blue squares) are plotted at the mean 
geometric radius for the sample, R = 0?20 (^ 60 kpc). 
For comparison, we also show measurements for the CC 
system in the inner regions of M87 from McLaughlin 
et al. (1994). For both the CCs and UCDs, the mea¬ 
surements are based on the method of moments formal¬ 
ism described in § 4.1 of McLaughlin et al. (1994). The 
CC system in M87, like the UCDs, is also quite similar 
to that of the galaxy in terms of flattening and position 
angle, a result which has been noted many times in the 
past (e.g., McLaughlin et al. 1994; Forte et al. 2012). 


3.3. Color Distributions 


• candidate globular cluster (g ^21.5) 

□ candidate UCD (u‘gz) + Virgo cluster member (B ^18) 



Fig. 20.— Same as Figure 18 except for the M60 region. The 
other galaxy whose GC system is visible at a ~ 190? 5 and (5 ~ 
11?6 is the massive, early-type galaxy M59. The spectroscopically 
confirmed UCDs in this field are the bright objects discussed in 
Strader et al. (2013), Seth et al. (2014) and Chilingarian & Mamon 
(2008). 



-0.5 0 0.5 -0.5 0 0.5 

Aa (deg) Aa (deg) 


Fig. 21. — Spatial distribution of bright GC- (gray dots) and 
UCD-like objects (open squares) in the four NGVS background 
fields. The star symbols plotted in the BG4 field show probable 
members of the 2 ; ~ 0.02 galaxy group identified by Crook et al. 
(2007). They are not identified as GCs or UCDs in our analysis. 


For the comparatively rich UCD system surrounding 
M87, our multi-band imaging — including Kg data from 
the NGVS-IR — makes it possible to carry out a com¬ 
prehensive study of UCD colors. As is well known, the 
GC systems of most massive galaxies are known to show 
multimodal — usually bimodal — color distributions (e.g. 
Ashman & Zepf 1992; Geisler et al. 1996; Gebhardt & 
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Fig. 22.— Radial variation in ellipticity, e, and position angle, 

for UCDs surrounding M87 (blue squares). For comparison, we 
also show profiles for the M87 halo light (dashed curve) and the 
globular cluster system from McLaughlin et al. (1994) (red circles). 

Kissler-Patig 1999; Kundu & Whitmore 2001; Kundu & 
Zepf 2007; Larsen et al. 2001a,b; Harris et al. 2006; Peng 
et al. 2006, 2009; Liu et al. 2011). This bimodality has 
been interpreted as evidence for a variety of physical pro¬ 
cesses that may have been involved in the formation of 
the host galaxy, including hierarchical growth (e.g., Cote 
et al. 1998, 2000, 2002; Beasley et al. 2002; Tonini 2013) 
and cluster formation in gas-rich major mergers (Ash¬ 
man & Zepf 1992; Muratov & Gnedin 2010). Alterna¬ 
tively, it has been proposed that a non-linear relationship 
between metallicity and (some) broadband color indices 
(Yoon et al. 2006, 2011, 2013; Cantiello & Blakeslee 2007; 
Blakeslee et al. 2012; Chies-Santos et al. 2012) is chiefly 
responsible for the observed bimodality. Once again, one 
should keep in mind that these processes need not be 
mutually exclusive, and it is likely that a combination 
of effects conspire to produce the ubiquitous bimodality 
observed for GC systems. 

Broadly speaking, existing photometry suggests that 
UCDs can also have a broad distribution in color, of¬ 
ten spanning the full range in GC colors (e.g., Hasegan 
et al. 2005; Mieske et al. 2006b; Evstigneeva et al. 2008; 
Madrid et al. 2010; Madrid 2011; Chiboucas et al. 2011). 
There is also some evidence that, at the highest lumi¬ 
nosities, this distribution gives way to a trend between 
luminosity and color, in the sense that the brightest sys¬ 
tems have the reddest colors (e.g., Mieske et al. 2004a; 
Blakeslee & Barber DeGraaff 2008), although there are 
at least some counterexamples of bright systems that are 
blue in color (Penny et al. 2014). 

Figures 23 shows optical/IR color histograms for our 
cleaned sample of 92 likely or confirmed UCDs in the 
M87 region. The UCD color distributions, plotted as the 
black histogram in each panel, include many of the pos¬ 
sible color combinations, from {u — g)o to {z — Ks)o. It 
is obvious from this figure that nearly all of the distri¬ 
butions are bimodal. To quantify this visual impression, 
we have used the Gaussian Mixture Modeling (GMM) 
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Fig. 23. — Optical and IR {u""grizKg) color distributions for 
our sample of 92 probable or confirmed UCDs in the M87 region. 
The blue and red curves show the individual components of the 
double Gaussian that best fits each distribution, determined using 
the GMM code assuming heteroscedasticity (see also Table 9). The 
dashed vertical lines show the inferred blue and red peaks, while the 
smooth black curve shows the sum of the individual components. 
The green histograms present the color distribution for the full 
sample of 64 dwarf nuclei. 


code of Muratov & Gnedin (2010) to evaluate the sig¬ 
nificance and characteristics of this apparent bimodality. 
Note that GMM is similar to the KMM code of Ashman 
et al. (1994), but it implements maximum likelihood pa¬ 
rameter estimation. The best-fit GMM models (blue and 
red curves in Figure 23) were calculated assuming het¬ 
eroscedasticity, i.e., independent dispersions for the two 
subcomponents. For comparison, the green histogram in 
each of these two figures shows the color distribution for 
the sample of 64 dwarf nuclei in this region. While a 
detailed analysis of the nuclei colors will be presented in 
a future paper in this series, we note that the nuclei, like 
the UGDs, are primarily blue. Unlike the UCDs, though, 
they show no strong evidence for bimodality. 

The best-fit, two-component Gaussians from our GMM 
analysis are shown in Figures 23 as the blue and red 
curves in each panel; their sum is shown in black. The 
two vertical dashed lines in each panel show the locations 
of the calculated blue and red peaks. The GMM-fitting 
results for the UCDs are summarized in Table 9. From 
left to right, this table records for each index, the mean 
color, /ia, and dispersion, da, as well as the mean color 
and dispersion for the separate blue (/i^, a^) and red (/i^, 
(Jr) components. The number of blue and red UCDs, 
Ni) and W, along with the overall fraction of blue ob¬ 
jects, / 5 , are given in the following three columns. The 
errors on each parameter were determined using a non- 
parametric bootstrap procedure. The final two columns 
give the D parameter and p-value for chi‘^ (see Muratov 
& Gnedin 2010 for details). As described in Muratov & 
Gnedin 2010, the D parameter should be larger than 2 
and p-value should be smaller than 0.05 for a bimodal 
distribution. As apparent from Figure 23, the color dis¬ 
tributions are found to be bimodal in all indices with 
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Fig. 24. — The measured half-light radius (r^) versus (g — Kg) 
color for our probable or coufirmed UCDs. Our selectiou criteria 
require that a UCD caudidate have r/j, > 11 pc. As showu iu Fig¬ 
ure 23, the vast majority of UCDs are blue iu color. The red UCDs 
are more compact, aud all have < 18 pc. It is possible that these 
red UCDs are simply iu the tail of the red GC size distributiou. 
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Fig. 25.— Optical giz) color distributious for UCD caudidates 
(black histogram) aud dwarf uuclei (greeu histogram) iu the M49 
regiou. 


the exceptions of (r — i)o and (i — z)o, both of which 
have p-values larger than 0.05. The lack of bimodality in 
these indices is perhaps not surprising given their limited 
wavelength baseline. 

The fraction of blue UCDs is found to be remarkably 
constant across color index, with a mean value of = 
0.85 ± 0.14. These estimates compare reasonably well 
with the value of 0.71 measured for GCs using the 
(g—z)o index for the 11 ACSVCS galaxies that fall within 
the M87 region considered here (Peng et al. 2006; Jordan 


et al. 2009). In other words, the UCD system in the core 
of Virgo (sub-cluster A) has a similar, or perhaps slightly 
larger, fraction of its members belonging to the blue sub¬ 
population as does the GC system in the same region. 

The interpretation of the UCD color distribution, how¬ 
ever, is even less straightforward than it is for GCs. Fig¬ 
ure 24 plots Th versus {g — Kg) to show that while the 
blue UCDs span a large range in sizes—from our lower 
limit of 11 pc to a maximum of ^ 40 pc—the red UCDs 
are all relatively compact. None of the red UCDs have 
larger than 18 pc. This raises the possibility that these 
red UCDs are simply the tail of the red GC population 
to large sizes. This does not mean, however, that there 
are no red UCDs. Some of the most massive confirmed 
UCDs in Virgo, such as M60-UCD1 (Strader et al. 2013; 
Seth et al. 2014) and M59c0 (Chilingarian & Mamon 
2008), have colors as red as the reddest GCs. Because 
of this ambiguity, it is difficult to interpret the color dis¬ 
tributions of UCDs. These issues will be explored for a 
variety of stellar systems (GCs, UCDs, galaxies, nuclei) 
in a future NGVS paper. For the time being, we simply 
note that UCDs as a population, defined using our selec¬ 
tion criteria, show bimodal color distributions, whereas 
dwarf nuclei are predominantly blue. The size distri¬ 
butions of the color subpopulations is consistent with a 
picture where some of the blue UCDs may have origins 
as stellar nuclei, and the red UCDs in our sample are 
possibly the largest red GCs. 

Figure 25 shows the optical color distribution of UCD 
candidates (black lines) and dwarf nuclei (green lines) 
in the M49 region. The six color indexes shown in this 
figure are based on the u^giz wavebands that we have 
in this region. The UCD colors do not show significant 
bimodal distributions, but most have “tails” on the red 
side, with the exception of the {g — i)o color index. Due 
to the limited number of UCD candidates in the M49 
region, we did not fit the color distribution using GMM 
pipeline. The green histograms show the distribution in 
color of the dwarf nuclei. The nuclei color distribution 
in the M49 region is similar to that in the M87 region: 
i.e., similar in color range with the UCDs and showing 
no obvious bimodality. 

3.4. Color Gradients 

As is well known, massive galaxies usually show ra¬ 
dial color gradients, in the sense that their colors become 
bluer with increasing distance from the center (i.e., a neg¬ 
ative gradient). This behavior is found in both early- and 
late-type galaxies (e.g., Boroson et al. 1983; Franx et al. 
1989; Peletier et al. 1990; de Jong 1996; MacArthur et al. 
2004; Liu et al. 2009; Gonzalez-Perez et al. 2011; Peletier 
et al. 2012). It is also known that the GC systems associ¬ 
ated with massive galaxies show negative color gradients 
(e.g., Strom et al. 1981; Geisler et al. 1996; Harris 2009a; 
Liu et al. 2011). In this case, the fairly strong color gra¬ 
dients are known to be driven mainly by the changing 
ratio of blue-to-red GCs as a function of radius, with the 
blue clusters being more spatially extended than their red 
counterparts (Geisler et al. 1996; Lee et al. 1998; Kundu 
et al. 1999; Larsen et al. 2001a; Peng et al. 2006, 2008). 

The dependence of UGD color on projected distance 
from M87 is shown in the upper panel of Figure 26. The 
sample consists of the 92 likely or confirmed UGDs (af¬ 
ter the rejection of spectroscopic non-members). UCDs 























































18 


Liu et al. 


1.8 

1.6 
o 1.4 

? 1.2 
bX) 

"" 1.0 
0.8 
0.6 

1.6 

1.4 

^ 1.2 
N 

U 1.0 

"" 0.8 
0.6 
0.4 

0.1 1.0 10.0 100.0 
R/Re 

Fig. 26.— Radial variation in color for UCDs and GCs in the M87 
region (upper and lower panels, respectively). The GCs show a 
significant negative color gradient (dashed line in the lower panel), 
becoming systematically bluer with increasing radius due mainly 
to the changing ratio of red-to-blue GCs. The blue and red GCs 
individually show shallower gradients (dotted lines in the upper and 
lower panels). The UCDs, on the other hand, show no evidence 
for a gradient in color (dashed line in the upper panel). The red 
circles in the upper panel show the dwarf nuclei in this region. 
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by the filled and open squares, respectively. The dashed 
line shows the best fit linear relation: 

(g - z)o = 0.876(±0.026) - 0.003(±0.024) log R/Re 

( 4 ) 

For comparison, we also show the nucleated dwarf galax¬ 
ies in the M87 region. Neither the dwarf nuclei, which 
he in the range 10 < R/Re ^ 50, nor the UCDs exhibit 
any trend with radius. 

The lower panel of Figure 26 shows the same relation 
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27.— Same as figire 26, but for the M49 region. 

and without, measured radial velocities are shown 
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for the GCs of M87. The dashed line in this case 

{g - z)o = 1.056(±0.003) - 0.193(±0.004) log R/Re 

( 5 ) 

refers to the entire GC system, with no division on color. 
Gonsistent with earlier findings (e.g. Jordan et al. 2004; 
Gantiello et al. 2007), we find that the varying ratio of 
blue to red GGs is chiefly responsible for the strong neg¬ 
ative gradient. The weaker trends observed for the sep¬ 
arate red and blue GCs, divided at {g — z)o = 1.00, are 
shown by the dotted lines: 

{g - z)red = 1.215(±0.007) - 0.068(±0.007) log R/Re 
{g-z)uue = 0.858(±0.003)-0.046(±0.003)logi?/i?e 

(6) 

For comparison, we overlay the fitted relations for the 
blue and red GCs on the UCDs in the upper panel of 
this figure. 

The existence of a blue tilt may bias the gradient mea¬ 
sured for the luminous blue GCs. To avoid this possi¬ 
bility, it is better to compare the UCDs and blue GCs 
in the same color and magnitude ranges. We therefore 
focus on the UCDs and GCs brighter than = 21.5 and 
bluer than {g — z)o = 1.0. The best fitted gradients in 
this case are V(^ — 2:)biueGC = —0.050 ±0.008 dex“^ and 
V(^ — ^)biueUCD = —0.018 ±0.018 dex“^. It is clear that 
luminous blue GC system has a significant gradient (at 
the ^ 6cr level) while blue UCD system shows no gradient 
at all. 

Figure 27 shows color gradients for GCs, UCDs and 
dwarf nuclei in the M49 region. The GCs in this figure 
are again only the luminous GCs which are brighter than 
g = 21.5. The best fitted lines (dashed lines) of the whole 
UCD system and whole GC system are: 

{g - z)ucD = 0.870(±0.050) - 0.004(±0.046) log R/Re 
{g - z)gc = 0.970(±0.020) - 0.082(±0.019) log R/Re 

( 7 ) 

If we fit only the UCDs and GCs that are brighter than 
g = 21.5 and bluer than {g — z)o = 1.0 , the gradi¬ 
ents are \/{g — 2 :)biueGC,M 49 = —0.038 ± 0.012 dex“^ and 
V(5'— 2 :)biueUGD,M 49 = —0.026 ±0.040 dex“^. Again, the 
luminous blue GC system shows a significant gradient 
(^3cr) while the blue UCD system shows no evidence of 
gradient. 

From the above analysis, it is clear that the GC sys¬ 
tems around M87 and M49 show significant color gra¬ 
dients, while the UCD systems show no gradient. We 
also notice that the gradients of UCDs have large errors, 
which makes the gradient values of UCDs and GCs com¬ 
patible at the ^ 2(7 level. A possible difference in the 
color gradients of GCs and UCDs may point to different 
formation mechanisms for the two classes of objects. 

3.5. Color Magnitude Relations 

The color-magnitude diagram for sources in the M87 
and M49 regions are shown in Figure 28 and 29. GCs are 
shown as gray dots, spectroscopically confirmed UCDs 
as filled squares, and candidate UCDs as open squares. 
Note that the GC sample has been cleaned of any objects 
that we have identified here to be UCDs on the basis of 
size. For reference, the 64 nuclei of dwarf galaxies in the 
M87 region and 22 nuclei in the M49 region are shown 
as red circles (open and filled). 
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Fig. 28.— Color magnitude diagram for the M87 region showing 
GCs, UCDs and dwarf nuclei. The blue and red sequences show 
the variations in the mean color for the blue and red GC sub¬ 
populations (see text for details). 



Fig. 29.— Color magnitude diagram for the M49 region. The 
UCDs (squares), dwarf nuclei (circles) and n^^iz-selected bright 
{go < 21.5) GC candidates (dots) are shown in this figure. 


As noted previously, the GC system of M87 (Mieske 
et al. 2006b; Forte et al. 2007; Peng et al. 2009; Harris 
2009b) exhibits a “blue tilt” (Harris et al. 2006; Strader 
et al. 2006; Mieske et al. 2006b), with the GCs belonging 
to the blue sub-population becoming progressively bluer 
with decreasing luminosity. To quantify this behavior, 
we bin the GCs by magnitude, and fit two-component 
Gaussians to the data points in each bin, thereby en¬ 
suring that all points are independent. The results are 
shown by the thin blue and red lines, which trace out 
the central locations of the blue and red components as 
a function of zo-band magnitude. Our analysis — which 
is based on a GC sample that is far larger than that used 
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Fig. 30.— The luminosity function of globular clusters in the M87 
region from the ACSVCS (black squares), -selected globular 

clusters from the NGVS {g < 24 mag, red circles) and the combined 
sample of NGVS globular clusters and UCDs (blue triangles). The 
dashed curve shows the gaussian luminosity function fit of Villegas 
et al. (2010) for M87. 

in any previous study of this galaxy — confirms the exis¬ 
tence of a significant blue tilt in the M87 GC system. For 
the M49 region, we only present bright GCs {g < 21.5) 
in the Figure 29. The blue line is the best fit for the blue 
GCs with [g — z){) < 1.0. Similar with M87 GC system, 
luminous blue GCs in the M49 region show significant 
blue tilt as well. 

Both Figure 28 and Figure 29 show that the UCDs 
participate in the general trends defined by the GCs. In¬ 
deed, at bright magnitudes, a significant fraction of the 
sources are actually UCD candidates, and so they have 
been at least partly responsible for driving the color- 
magnitude trends noted in previous ground-based stud¬ 
ies. The abrupt bend to redder colors exhibited by the 
UCDs is similar to the behavior of UCD candidates in 
Abell 1689 noted by Mieske et al. (2004b), although 
UCD membership constraints in that case were fairly 
limited. For both M87 and M49 regions, many of the 
bright UCDs are confirmed spectroscopic members (filled 
squares), leaving little doubt over the reality of this color- 
magnitude relation. Interestingly, the dwarf nuclei seem 
to follow the trends defined by the UCDs and (blue) GCs. 
It also noteworthy that all of the 11 UCDs in the M87 re¬ 
gion and all of the 3 UCDs in the M49 region with colors 
redder than {g — z)o > 1.0 are brighter than zq — 20. 

3.6. Luminosity Functions 

The M87 GC luminosity function has been studied 
many times in the past (e.g., McLaughlin et al. 1994; 
Kundu et al. 1999; Jordan et al. 2002, 2006, 2007; Peng 
et al. 2009; Villegas et al. 2010), so it is possible to com¬ 
pare these earlier results (many of which are based on 
high-resolution imaging from HST) with our new mea¬ 
surements from the NGVS. Of particular interest in this 
context is a comparison of the UCDs to GCs to see if the 
two populations are consistent with being drawn from the 
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same parent distributions in luminosity. Although pre¬ 
vious studies of this sort have generally failed to provide 
conclusive evidence for, or against, the hypothesis that 
the UCDs and GCs are drawn from the same luminos¬ 
ity function (due largely to limitations in UCD sample 
size and the heterogeneity of the GC and UGD samples 
used for the comparisons; e.g., Mieske et al. 2004a, 2007; 
Gregg et al. 2009), some recent investigations have ar¬ 
gued that the UGD luminosities are consistent with their 
being the most luminous members of the host galaxy’s 
GG system (e.g., Hilker 2009; Norris & Kannappan 2011; 
Mieske et al. 2012). 

Figure 30 shows the GG luminosity function of M87 
according to the AGSVGS (black squares). This HST- 
based measurement has the advantage that typical GGs 
at the distance of Virgo can be resolved by HST/AGS 
(e.g., Jordan et al. 2005), thus significantly reducing con¬ 
tamination. The data have been taken from the catalog 
of Jordan et al. (2009), and consist of 1748 GGs identi¬ 
fied in a single AGS field located at the center of M87. 
To allow a direct comparison to the NGVS results, the 
^-band magnitudes from the AGSVGS have been trans¬ 
formed from the SDSS to MegaGam photometric system 
— a typical correction of ^ 0.078 mag in the sense that 
the MegaGam magnitudes are brighter^^. 

The dashed black curve in this figure shows a Gaus¬ 
sian fitted to the transformed AGSVGS data (Villegas 
et al. 2010). The red circles show the corresponding 
GG luminosity function from the NGVS. The luminosity 
functions were scaled vertically to match at the peak of 
the distribution (at g ^ 23) where the statistical errors 
are the smallest and both samples should be compete. 
Any objects in the NGVS GG catalog that have mea¬ 
sured radial velocities inconsistent with their being GGs 
{vr > 3500 km s“^) have also been removed, leaving a 
sample of 9497 objects. Finally, we add to the NGVS 
GG sample the 92 confirmed or probable UGDs from 
this study, and plot the combined luminosity function 
of these 9589 objects as blue open triangles. 

There is little to differentiate the two NGVS luminosity 
functions, although the latter sample (which includes the 
UGDs) is more noticeably enhanced over the Gaussian 
curve for ^0 ^ 21. It is difficult to say with certainty 
whether either NGVS sample provides a superior match 
to the AGSVGS results because the UGDs represent only 
a tiny enhancement (~ 1-2%) to the GG sample, making 
it is difficult to draw firm conclusions. These difficulties 
are exacerbated by the differences in the spatial coverage 
of the HST and NGVS surveys, combined with possible 
differences in the spatial distribution of GGs and UGDs. 

While our ability to isolate large and nearly 
contaminant-free sample of UGDs and GGs has improved 
considerably in recent years, thanks to both spectro¬ 
scopic surveys and imaging programs like the AGSVGS 
and NGVS, we continue to lack a fully developed and 
quantitative theory for the GG luminosity function (i.e., 
predictive models). Until then, definite conclusions on 
the correspondence between the UGD and GG luminos¬ 
ity functions will likely remain elusive. 

3.7. Connections to Features in the Intracluster Light? 

For more information, see http://www2.cadc-ccda.hia- 
iha.nrc-cnrc.gc.ca/megapipe/docs/filters.html 


The innermost regions of the principal sub-clusters in 
Virgo are crowded and dynamic environments in which 
galaxy interactions have likely played an important role 
in shaping their present-day appearance. Some UGD for¬ 
mation scenarios invoke processes like satellite stripping, 
accretion and disruption, so it is natural to ask if there is 
evidence for a correspondence between the UGDs found 
in this study and any previously identified features in the 
diffuse IGL. Given the modest size of the UGD system in 
M60, we confine the discussion to the UGDs surrounding 
M87 and M49. 

There is evidence for past interactions in both of these 
regions from several independent studies, including on¬ 
going dwarf galaxy accretion (McNamara et al. 1994; 
Lee et al. 1997; Arrigoni Battaia et al. 2012), possible 
kinematic substructuring among GGs and PNe (Gote 
et al. 2003; Doherty et al. 2009; Romanowsky et al. 2012; 
Agnello et al. 2014), dynamical structures (e.g., Boselli 
et al. 2014), indirect evidence for tidal stripping of cer¬ 
tain types of galaxies, including compact ellipticals (e.g., 
Gote et al. 2008; Gote 2010; Peng et al. 2008; Goenda 
& Muriel 2009), and the discovery of discrete features in 
the surrounding IGL (e.g., shells, streams and plumes) 
that almost certainly arose from galaxy interactions and 
accretions (Mihos et al. 2005, 2013; Rudick et al. 2009; 
Janowiecki et al. 2010; Ferrarese et al. 2012). 

Figure 31 shows some of these features for the M87 
(upper panels) and M49 (lower panels) regions. The 
left panels show the original U-band images from the 
study of Mihos et al. (2005), while the images on the 
right show residuals after subtracting the best-fit galaxy 
model. The irregular polygons in each panel indicate the 
features identified by Janowiecki et al. (2010). The blue 
and red crosses show the UGD candidates in each re¬ 
gion after dividing the sample by color at {g — z)q = 1.0. 
The UGD samples shown here are the cleaned sample 
of UGDs in M87 (see §2.6), and the selected UCD 
candidates in M49 (see Table 8). 

Of the five IGL features in M87 noted by Janowiecki 
et al. (2010), only eight UCDs fall within their bound¬ 
aries (a single object in region and seven more in 
region #4). In the case of M49, two regions contain one 
or more UCDs: two objects in region ^^1, and one more 
in region 7^5. All in all, though, there is little evidence 
that the UCD populations are preferentially associated 
with the catalogued features in the IGL of these sub¬ 
clusters. Of course, this may simply refiect the fact that 
most of the known IGL features are found at large radii 
(where there are relatively few UCDs) presumably be¬ 
cause phase-mixing timescales in the inner regions would 
rapidly erase any such features. 

4. DISCUSSION 

We showed in §2 that it is possible to measure reli¬ 
able sizes for UCDs larger than rh ^ 10 pc using NGVS 
images. By combining these size estimates with cuts in 
the color-color diagram, visual classifications, and new 
or published spectroscopic redshifts, we have produced 
relatively clean samples of UCDs for ^4 deg^ regions 
centered on each of M87, M49 and M60 — the three 
most massive galaxies in the Virgo cluster. In this sec¬ 
tion, we summarize our key findings and discuss their 
implications for UCD formation models. 

Perhaps most fundamentally, there is clear evidence for 
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Fig. 31. — (Upper panels). Mosaic U-band image of M87 from the study of Mihos et al. (2005). The red and blue UCD candidates 
identified from the NGVS are shown as the red and blue crosses. The cyan cross shows the M87 center. Regions outlined in magenta are 
the low surface brightness features identified by Janowiecki et al. (2010). The left and right panels show the original and model galaxy 
subtracted image, respectively. (Lower panels). The same as above, except for M49, and showing UCD candidates identified in this region. 


large region-to-region differences in the number of UCDs 
(§3.1). Such variations would not be expected if A/'ucd 
scales with either the luminosity or stellar mass of the 
central galaxy, which are similar in all three cases. How¬ 
ever, these variations would be expected if the number of 
UCDs scales in proportion to the mass of the surrounding 
subcluster; as Figure 17 shows, this is precisely what is 
observed, with N\jcb increasing in step with other trac¬ 
ers of sub-cluster mass: X-ray gas mass, dark matter 
mass, and the number of surrounding galaxies and GCs. 
Moreover, the flattening and orientation of the UCDs in 
M87 appear to match that of the central galaxy isophotes 
and/or ICL in this region. 

An increasing number of studies support the view that 
UCD formation is tied to tidal stripping (e.g., Bekki et al. 
2001; Drinkwater et al. 2003; Pfeffer & Baumgardt 2013; 
Pfeffer et al. 2014; Seth et al. 2014; Janz et al. 2015; 
Mihos et al. 2015). In this scenario, the larger UCD sys¬ 
tem surrounding M87 3.5 and 7.8 times more than 

M49 and MOO) can be understood in terms of the larger 
reservoir of low-mass galaxies in this more massive en¬ 


vironment, while the stronger tidal forces (due to the 
larger virial mass) would allow low-mass galaxies to be 
stripped more easily. 

Meanwhile, the rather blue colors of most UCDs (§3.3) 
are consistent with expectations in the tidal stripping 
scenario if their putative progenitors were low-mass, 
metal-poor galaxies. Similarly, the fact that both the 
UCD colors and their location in the color-magnitude re¬ 
lation closely match the nuclei lends further support to 
the notion that tidal stripping of low-mass galaxies was a 
dominant process in UCD formation. Although we find 
no evidence that the UCDs are preferentially associated 
with known low-surface brightness features in the ICL, 
this probably reflects the fact that only those features 
located at large radii — beyond the bulk of the UCD 
populations in each subcluster — would still be visible at 
the present time. Clearly, though, this is an area where 
additional study is needed, preferably using numerical 
simulations with realistic cosmological initial conditions. 

In a companion paper, Zhang et al. (2015) have exam¬ 
ined the kinematics of UCDs in the M87 region, com- 
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Fig. 32. — Mosaic of images showing 20 of the objectively- 
selected UCD candidates from Tables 2 and 3. The first nine panels 
show known nucleated dwarf galaxies (class=2) from Binggeli et al. 
(1985) or new discoveries from the NGVS. The objects have been 
loosely arranged according to the visual prominence of an under¬ 
lying diffuse halo, which is unmistakable in the first objects and 
non-existent in last ones. Eighteen of these 20 objects are con¬ 
firmed radial velocity members of the Virgo cluster, the exceptions 
being two dwarf nuclei (object IDs = 7 and 12). The measured 
velocities and the envelope parameters Aenv (see equation 8) are 
labeled at the bottom of each panel. 

paring their properties to GCs and low-mass galaxies. 
They find the UCDs have a similar velocity dispersion 
profile with the blue GCs. There is a hint that the ve¬ 
locity anisotropy of the UCDs changes from tangentially 
biased < 0) to radially biased (/3 > 0) with increas¬ 
ing distance from M87, and, to within the uncertainties, 
the GCs and UCDs are found to have similar anisotropy 
profiles with possible deviations at large radii. On the 
other hand, the UCDs show somewhat larger rotation 
than either the blue or red GCs. Thus, while there is 
no compelling evidence from kinematics alone that tidal 
stripping of low-mass galaxies is a primary mechanism for 
UCD formation, the more radially biased orbits of UCDs 
compared to blue GCs at large radii are suggestive that 
stripping plays a role. 

Is there evidence from the structural properties of the 
UCDs and dwarf galaxies themselves that might support 
a connection between UCDs and low-mass galaxies? Fig¬ 
ure 32 uses our sample of 92 likely or confirmed UCDs 
surrounding M87, along with the nine dwarf nuclei in 
this same region that were objectively identified in our 
UCD search, to illustrate a heuristic structural sequence 
among these ^100 objects. The 20 objects shown in this 
mosaic have been arranged, from left to right and top 
to bottom, by the apparent prominence of a diffuse en- 
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Fig. 33.— The surface brightness profiles (blue solid lines), fitted 
King models (cyan dashed lines) and PSFs (red dotted lines) of 20 
objectively-selected UCD candidates which are shown in Figure 32. 

velope surrounding each compact object. For the first 
half of the sample, there is clear evidence for the two- 
component structure (i.e., a nucleus plus an envelope) 
that is characteristic of dE,N, dS0,N, or dSph,N galax¬ 
ies. Indeed, all but one of these objects was previously 
cataloged as a nucleated dwarf galaxy residing in the 
cluster core (Binggeli et al. 1985). The sequence extends 
to include very compact UCD candidates that show lit¬ 
tle or no sign of an extended envelope. Figure 33 shows 
surface brightness profiles, PSFs and fitted King models 
for each of the objects shown in Figure 32. The gradual 
progression from a two- to one-component structure is 
apparent. 

In Figure 34, we plot the sequence of UCD and enve¬ 
lope structure on a common axis to highlight the change 
in envelope fraction relative to the central component. 
This figure is reminiscent of Figure 5 in Pfeffer & Baum- 
gardt (2013), which showed the simulated surface bright¬ 
ness profile evolution of a dwarf galaxy being tidally 
stripped. 

Because this sequence is based purely on visual ap¬ 
pearance, a more quantitative approach is desirable. We 
therefore introduce the parameter 

Aenv = ^16 — ^32 (8) 

where qiq and g ^2 are the ^-band corrected point-source 
magnitudes measured in apertures of diameter 16 and 
32 pixels (see §2.2). These apertures correspond to radii 
of 1?5 and 3?0, respectively, and their difference should 
provide a diagnostic of any diffuse light surrounding the 
candidates. 

Figure 35 shows the distribution of Aenv values for 
the sample of 92 likely or confirmed UCDs, and 64 
dwarf nuclei, in the M87 region. The dwarf nuclei. 
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Fig. 34.— Surface brightness profiles of UCDs with a range 
of envelope fraction. The surface brightness profiles were shifted 
vertically to match the profile of ID=48 at D ~ 30pc (roughly 
the mean half-light radius of UCDs) to highlight the differences in 
envelope components. 



Fig. 35.— Histogram of difference in aperture magnitudes, Aenv, 
for UCD candidates and dwarf nuclei in the M87 region. Here 
Aenv is the magnitude difference measured in apertures of 16- and 
32-pixel diameter. A value of Aenv = 0.05 mag (denoted by the 
dotted vertical line) roughly divides objects with, and without, 
diffuse extended envelopes. 

which are shown in red, are broadly distributed with 
0.00 < Aenv ^ 1-00 mag. By contrast, the UCDs (class 
= 1 objects), which are shown by the cyan and blue his¬ 
togram, are strongly peaked at Aenv — 0 mag, with an 
asymmetry to positive values. We can use the observed 
dwarf distribution as a guide in separating UCD can¬ 
didates with, and without, an appreciable envelope by 
adopting Aenv — 0.05 mag as the dividing point. UCD 


Fig. 36.— Cumulative distributions for Dmst (fhe projected dis¬ 
tance from M87) for the objectively identified UCD candidates in 
the M87 region (blue and cyan lines), the full sample of 64 dwarf 
nuclei (red solid line), and the nine additional candidates subse¬ 
quently classified as dwarf nuclei (red dotted line). The cleaned 
sample of 92 UCD candidates has been divided into two classes by 
Aenv, which is a measure of the prominence of an outer envelope. 
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Fig. 37.— Envelope prominence as a function of the projected 
distance from M87 for UCDs and nucleated dwarf galaxies. 

candidates with Aenv <0.05 mag and Aenv >0.05 mag 
are shown in blue and cyan, respectively. 

The cumulative distribution of projected distances to 
M87, Dm 875 foi* these three samples are compared in Fig¬ 
ure 36. It is obvious that there are dramatic differences 
in their spatial distribution, with the dwarf nuclei lying 
preferentially at larger distances than the UCD candi¬ 
dates (see Drinkwater et al. 2004 for a similar result in 
the Fornax cluster). UCDs that show signs of a faint en¬ 
velope, in turn, have a larger mean distance than those 
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4nv = gi6 - g32 (mag) 

Fig. 38.— The same as Figure 35, except for the M49 region. 



that show no envelope. These same trends are also ap¬ 
parent in Figure 37, which plots Aenv against Dmst foi* 
the individual objects. Although the dwarfs have a distri¬ 
bution in Aenv 5 most are located at large distances (Dmst 
> 0?2). The UCD candidates with Aenv < 0-05 mag are 
found at all radii, whereas those with Aenv ^ 0-05 mag 
appear to form an intermediate population. 

As a check on the significance of the observed trend for 
compact systems surrounding M87, we show in Figure 38, 
39 and 40 the analogous distribution for M49, but for 
its sample of 28 UCDs and 22 nucleated dwarf galaxies. 
Clearly, the behavior seen in the M87 region is seen here 
as well. 

The most straightforward interpretation of the ob- 
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Fig. 40.— The same as Figure 37, except for the M49 region. 

served trend is that proximity to the central, massive 
galaxy in these regions governs the envelope structure 
of UCDs and dwarf nuclei. This finding may provide 
an important new piece of evidence in support of the 
tidal stripping scenario, since the strong tidal forces near 
the bottom of the gravitational potential wells in these 
sub-clusters would act to strip away the envelopes of the 
nearest UCD progenitors. 

5. CONCLUSIONS AND DIRECTIONS FOR FUTURE WORK 

In this paper, we have presented the first homoge¬ 
neous and complete photometric study of UCDs around 
the three prominent Virgo sub-clusters, centered on the 
galaxies M87, M49, and M60. Using deep, multi-color 
NGVS imaging, we select UCDs with high completeness 
using their measured colors and half-light radii. Below, 
we summarize the main findings of our study: 

• Using the excellent image quality in the NGVS 
data, we can reliably measure the half-light radii of 
objects in Virgo down to a limiting size of ^ 10 pc 
at (7 ^ 21.5 mag, making the data ideally suited for 
UGD studies. 

• We present samples of 92, 28, and 23 confirmed 
or probable UGDs associated with M87, M49, and 
M60 (within 2° x 2° boxes), respectively, selected 
to have 18.5 < g < 21.5 mag (—12.7 < Mg < —9.7) 
and half-light radii 11 < rh < 100 pc. These ob¬ 
jects were selected using a combination of criteria 
including luminosity, size, colors, effective surface 
brightness, and visual inspection. 

• The number of UCDs scales with the total mass of 
the host system. M87 contains ^ 3.5 and 7.8 times 
more UGDs than M49 and MOO within a radius of 
25i?e, respectively. However, there are tight corre¬ 
lations between the number of UGDs and the pa¬ 
rameters that describe the total mass of each sub¬ 
cluster (X-ray gas mass, dark matter mass, num¬ 
ber of neighboring galaxies, and number of GGs). 
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Normalized to the parameters such as dark matter 
mass and the number of neighboring galaxies, M87, 
M49 and M60 appear to have formed UCDs with 
comparable efficiency. 

The spatial distribution of UCDs mirrors the galaxy 
light. The M87 UCDs appear to follow the ellip- 
ticity and orientation of the low surface brightness 
stellar light in the halo. 

Most UCDs have blue eolors. We find that ^85% 
of the UCDs in M87 have colors similar to blue 
CCs or stellar nuclei. The remaining red UCDs 
are mostly more compact (11 < r^ < 18 pc, and 
may be consistent with being the tail of the red 
CC population. However, we note that a few of the 
most luminous UCDs in the literature are quite red 
in color, so not all UCDs are blue. 

M87 and M49 UCDs show no eolor gradient with 
galaetoeentrie radius. The UCDs in the M87 and 
M49 regions display no significant gradient in color. 
This lack of a gradient is also found for the nuclei of 
dwarf galaxies around M87 and M49. Meanwhile, 
the CCs show a color gradient in both their blue 
and red sub-populations. 

M87 and M49 UCDs show a relation between eolor 
and magnitude. Like what has been seen 

before for blue CCs, the M87 UCDs display a 
color-magnitude relation, where the more luminous 
UCDs tend to be redder. Dwarf galaxy nuclei also 
seem to follow a similar color-magnitude trend. 
Similar results are found in the M49 region. 

The M87 UCD system represents a 1-2% enhanee- 
ment over a Gaussian GCLF at the bright end. By 
combining the GCLF from the ACS Virgo Cluster 
Survey and the NGVS, we find that inclusion of the 
UCDs causes the luminosity function at the bright 
end to be enhanced over a Gaussian extrapolation 
by a small amount. Lacking a developed theory 
of the GCLF, however, the interpretation of this 
result is difficult. 

There is little evidenee that UCDs are preferentially 
assoeiated with streams or features eurrently visi¬ 
ble in the ICL. We matched the locations of UCDs 
with ICL features detected in deep images of the 
M87 and M49 regions and did not find that UCD 
numbers were enhanced in these regions. 

UCDs around M87 follow a morphologieal sequenee 
defined by the prominenee of their outer, low sur- 
faee brightness envelope, with the envelope frae- 
tion eorrelating with distanee from M87. We quan¬ 
tify the fraction of stellar light in UCDs that is 
contained in an outer envelope, and find that the 
M87 UCDs span a range of envelope fraction that 
merges into the sequence of nucleated dwarf galax¬ 
ies. The mean UCD envelope fraction rises with 
distance from M87. The similar result is also found 
for the UCDs around M49. This sequence suggests 
that tidal stripping may play an important role in 
the formation of UCDs. 


Taken together, our results suggest that many UCDs 
are distinct from normal CCs, and originate as stellar nu¬ 
clei that are subsequently subjected to a process of tidal 
stripping that depends mainly on the total mass of the 
host halo. Some objects selected as UCDs in our study, 
however, may be the massive and extended tail of the 
GC population as these two populations likely overlap. 

Several obvious extensions to this work present them¬ 
selves. First, our study has focused on just three regions 
covering a combined area of ^ 11 deg^. Since the full 
NGVS survey spans an area of ^ 104 deg^ (Ferrarese 
et al. 2012), a systematic search for UCDs — based on 
both effective radii and location in the u'^gz diagram — 
covering a ten-fold larger area is possible and, indeed, will 
be presented in a future paper in this series (Liu et ah, 
in preparation). Second, the addition of i^g-band imag¬ 
ing is obviously invaluable in the identification of both 
UCDs and CCs (Munoz et al. 2014; §2.4.1) as it leads 
directly to a 2-3 x reduction in the level of background 
contamination (see §2.5). Programs to acquire additional 
IR imaging are currently underway (CFHT/WIRCam, 
VISTA/VIRCam). 

Third, a careful examination of the scaling relations 
for UCDs and CCs requires accurate structural infor¬ 
mation (e.g, effective radius, mean surface brightness, 
concentration, etc.) that is beyond the reach of ground- 
based imaging. While the NGVS image quality is ex¬ 
ceptional compared to most ground-based surveys, it is 
nevertheless suitable only for identifying compact stellar 
systems with sizes larger than rh ^ 10 pc (see §2.3). On 
the other hand, measuring accurate structural parame¬ 
ters for Virgo UCDs and CCs is entirely straightforward 
with HST, and both new and archival observations could 
be used for this purpose. Fourth, low-resolution spec¬ 
troscopy for UCD candidates from the NGVS is within 
the reach of existing multi-fiber spectrographs on 4m- 
and 8m-class telescopes, and such spectra would allow 
the kinematic and chemical properties of NGVS-selected 
UCDs to be examined in a systematic way (e.g., Zhang 
et al. 2015). Finally, high-resolution, integrated-light 
spectroscopy for the brightest UCDs would, once com¬ 
bined with structural parameters from HST and multi¬ 
color photometry from the NGVS and UV/IR surveys, 
lead to a better understanding of their internal proper¬ 
ties, including their initial mass function (Dabringhausen 
et al. 2010, 2012; Marks et al. 2012), dynamical evolu¬ 
tion (Chilingarian et al. 2011), and possible dark matter 
and/or central black hole content (Hasegan et al. 2005; 
Mieske et al. 2008, 2013; Frank et al. 2011; Chilingarian 
et al. 2011; Seth et al. 2014). 
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TABLE 4 

Photometric Properties of UCD Candidates in the M49 Region 


ID 

(1) 

NGVSID 

(2) 

«J2000 

(deg) 

(3) 

^J2000 

(deg) 

(4) 

E(B - V) 
(mag) 

(5) 

516 

(mag) 

(6) 

('«'* - 5)8 
(mag) 

(7) 

(5 - i)8 
(mag) 

(8) 

(9 - 
(mag) 

(9) 

^env 

(mag) 

(10) 

UGD 

(11) 

1 

NGVS-J122708.94 + 074228.3 

186.7872607 

7.7078639 

0.024 

19.429 + 0.001 

1.025 + 0.003 

0.813 + 0.002 

0.955 + 0.002 

0.239 

1 

2 

NG VS-J122718.40 + 081154.2 

186.8266527 

8.1983810 

0.023 

20.966 + 0.002 

1.076 + 0.008 

0.794 + 0.004 

0.900 + 0.007 

-0.140 

0 

3 

NG VS-J122735.11 + 072518.7 

186.8962982 

7.4218483 

0.022 

20.992 + 0.003 

0.923 + 0.008 

0.678 + 0.005 

0.790 + 0.008 

0.085 

0 

4 

NGVS-J122741.97+075613.3 

186.9248642 

7.9370367 

0.025 

21.193 + 0.003 

1.146 + 0.012 

0.747+0.006 

0.862 + 0.010 

0.075 

0 

5 

NG VS-J122744.26 + 072542.5 

186.9344272 

7.4284640 

0.022 

20.645 + 0.002 

0.968 + 0.006 

0.683 + 0.004 

0.778 + 0.006 

0.103 

1 

6 

NGVS-J122803.46 + 065343.1 

187.0144136 

6.8953085 

0.020 

21.446 + 0.004 

0.987+0.011 

0.737+0.007 

0.868 + 0.011 

0.140 

1 

7 

NGVS-J122806.14 + 065909.0 

187.0255686 

6.9858243 

0.020 

20.272 + 0.002 

1.485 + 0.008 

0.984 + 0.004 

1.250 + 0.005 

0.321 

1 

8 

NGVS-J122806.31 + 080342.8 

187.0262901 

8.0618891 

0.022 

21.475 + 0.003 

0.923 + 0.013 

0.697+0.007 

0.795 + 0.013 

0.092 

0 

9 

NGVS-J122818.37+074924.9 

187.0765475 

7.8235740 

0.024 

21.563 + 0.003 

1.011 + 0.013 

0.787+0.007 

0.935 + 0.012 

0.148 

0 

10 

NGVS-J122828.41 + 083115.5 

187.1183613 

8.5209855 

0.018 

20.961 + 0.002 

0.984 + 0.007 

0.707+0.005 

0.832 + 0.007 

0.124 

0 

11 

NGVS-J122828.60 + 070228.2 

187.1191854 

7.0411563 

0.018 

20.713 + 0.002 

0.907+0.007 

0.799 + 0.005 

0.892 + 0.007 

0.221 

1 

12 

NGVS-J122830.49 + 071929.1 

187.1270287 

7.3247459 

0.021 

20.459 + 0.002 

1.113 + 0.006 

0.795 + 0.004 

0.947 + 0.006 

0.425 

0 

13 

NG VS-J122849.25 + 075919.4 

187.2051905 

7.9887156 

0.022 

21.260 + 0.003 

0.918 + 0.010 

0.844 + 0.006 

0.914 + 0.010 

0.077 

1 

14 

NG VS-J122850.66 + 082621.0 

187.2110768 

8.4391793 

0.021 

21.534 + 0.004 

0.930 + 0.013 

0.759 + 0.007 

0.864 + 0.013 

0.069 

0 

15 

NGVS-J122857.99 + 083513.8 

187.2416238 

8.5871653 

0.021 

18.812 + 0.001 

1.479 + 0.003 

1.072 + 0.001 

1.325 + 0.001 

0.184 

0 

16 

NGVS-J122858.00 + 071759.9 

187.2416690 

7.2999606 

0.022 

20.185 + 0.001 

1.119 + 0.005 

0.892 + 0.003 

1.056 + 0.004 

0.154 

0 

17 

NGVS-J122906.90 + 073558.8 

187.2787298 

7.5996681 

0.023 

21.230 + 0.003 

0.986 + 0.009 

0.759 + 0.006 

0.827 + 0.010 

0.129 

0 

18 

NGVS-J122917.85 + 082558.1 

187.3243594 

8.4328050 

0.022 

18.640 + 0.001 

1.320 + 0.003 

0.993 + 0.001 

1.217 + 0.002 

0.185 

0 

19 

NGVS-J122922.85 + 074901.2 

187.3452098 

7.8169941 

0.025 

20.813 + 0.002 

1.022 + 0.007 

0.762 + 0.004 

0.843 + 0.006 

0.088 

1 

20 

NGVS-J122924.89 + 072105.3 

187.3536954 

7.3514762 

0.020 

21.317 + 0.003 

0.917 + 0.010 

0.670 + 0.008 

0.721 + 0.012 

0.244 

0 

21 

NGVS-J122925.10 + 075407.0 

187.3545827 

7.9019461 

0.024 

20.231 + 0.001 

1.030 + 0.005 

0.762 + 0.003 

0.842 + 0.004 

0.051 

1 

22 

NGVS-J122926.23 + 081658.8 

187.3593098 

8.2829991 

0.022 

21.453 + 0.003 

1.076 + 0.012 

0.746 + 0.007 

0.892 + 0.012 

0.036 

1 

23 

NG VS-J122928.21 + 075355.4 

187.3675253 

7.8987154 

0.024 

19.860 + 0.001 

1.266 + 0.005 

0.918 + 0.002 

1.030 + 0.003 

0.093 

1 

24 

NG VS-J122928.62 + 075824.6 

187.3692499 

7.9734990 

0.023 

20.932 + 0.004 

1.123 + 0.012 

0.720 + 0.006 

0.819 + 0.011 

0.033 

1 

25 

NG VS-J122934.41 + 080540.9 

187.3933684 

8.0947073 

0.022 

20.493 + 0.002 

0.977 + 0.006 

0.725 + 0.004 

0.804 + 0.006 

0.092 

0 

26 

NG VS-J122948.37+080042.0 

187.4515565 

8.0116645 

0.022 

20.741 + 0.004 

1.084 + 0.010 

0.808 + 0.006 

0.930 + 0.009 

0.105 

1 

27 

NGVS-J122957.08 + 074820.0 

187.4878332 

7.8055582 

0.024 

19.737 + 0.001 

1.035 + 0.004 

0.851 + 0.002 

0.985 + 0.003 

0.063 

1 

28 

NGVS-J122958.06 + 080409.2 

187.4919199 

8.0692203 

0.022 

19.754 + 0.001 

1.197 + 0.004 

0.896 + 0.002 

1.041 + 0.003 

0.016 

1 

29 

NGVS-J123002.49 + 075846.7 

187.5103664 

7.9796416 

0.022 

20.983 + 0.003 

1.004 + 0.010 

0.738 + 0.006 

0.833 + 0.010 

0.064 

1 

30 

NGVS-J123004.35 + 073932.2 

187.5181421 

7.6589449 

0.020 

21.511 + 0.004 

0.967 + 0.010 

0.716 + 0.008 

0.799 + 0.011 

0.079 

1 

31 

NGVS-J123006.03 + 074901.2 

187.5251438 

7.8169882 

0.023 

21.255 + 0.003 

1.003 + 0.009 

0.738 + 0.006 

0.879 + 0.009 

0.056 

1 

32 

NGVS-J123006.89 + 080509.5 

187.5287169 

8.0859729 

0.022 

20.455 + 0.002 

0.999 + 0.005 

0.769 + 0.003 

0.856 + 0.005 

0.069 

1 

33 

NGVS-J123008.39 + 082507.4 

187.5349653 

8.4187179 

0.021 

19.044 + 0.001 

1.210 + 0.002 

0.873 + 0.001 

1.044 + 0.002 

0.212 

1 

34 

NGVS-J123008.69 + 065353.4 

187.5362036 

6.8981683 

0.017 

20.164 + 0.001 

0.982 + 0.004 

0.818 + 0.003 

0.947 + 0.004 

0.109 

0 

35 

NGVS-J123009.17+074127.5 

187.5381980 

7.6909633 

0.021 

19.526 + 0.001 

1.316 + 0.003 

0.864 + 0.002 

1.067 + 0.002 

0.086 

1 

36 

NGVS-J123015.56 + 083445.0 

187.5648407 

8.5791729 

0.020 

20.896 + 0.002 

1.258 + 0.007 

0.887+0.004 

1.043 + 0.006 

0.049 

1 

37 

NGVS-J123016.65 + 080454.4 

187.5693858 

8.0817884 

0.022 

20.859 + 0.002 

1.035 + 0.007 

0.801 + 0.004 

0.891 + 0.006 

0.074 

1 

38 

NGVS-J123018.06 + 075024.9 

187.5752391 

7.8402478 

0.021 

19.466 + 0.001 

1.175 + 0.003 

0.862 + 0.002 

1.045 + 0.002 

0.210 

1 

39 

NGVS-J123029.94 + 074755.5 

187.6247355 

7.7987386 

0.021 

20.397 + 0.002 

1.185 + 0.005 

0.831 + 0.003 

0.996 + 0.004 

0.053 

1 

40 

NGVS-J123047.31 + 073619.4 

187.6971300 

7.6053799 

0.019 

20.161 + 0.001 

0.982 + 0.005 

0.752 + 0.003 

0.878 + 0.004 

0.282 

0 

41 

NGVS-J123048.49 + 065420.2 

187.7020501 

6.9056234 

0.017 

21.457 + 0.003 

1.144 + 0.011 

0.774 + 0.006 

0.905 + 0.010 

0.094 

0 

42 

NGVS-J123109.91 + 080405.3 

187.7912805 

8.0681353 

0.024 

19.812 + 0.001 

0.977 + 0.003 

0.744 + 0.002 

0.865 + 0.003 

0.124 

1 

43 

NGVS-J123113.76 + 074811.9 

187.8073500 

7.8033039 

0.020 

21.293 + 0.003 

0.993 + 0.009 

0.745 + 0.006 

0.843 + 0.009 

0.056 

1 

44 

NGVS-J123148.91 + 071620.9 

187.9537817 

7.2724780 

0.019 

19.244 + 0.001 

1.237 + 0.003 

0.803 + 0.002 

1.032 + 0.002 

0.205 

0 

45 

NGVS-J123152.37+064659.4 

187.9681966 

6.7831568 

0.021 

20.163 + 0.001 

0.972 + 0.005 

0.778 + 0.003 

0.900 + 0.004 

0.157 

0 

46 

NG VS-J123212.59 + 074047.6 

188.0524452 

7.6798850 

0.021 

21.226 + 0.003 

0.944 + 0.009 

0.716 + 0.006 

0.839 + 0.009 

0.098 

0 

47 

NG VS-J123320.46 + 080226.6 

188.3352586 

8.0407230 

0.019 

20.996 + 0.002 

1.011 + 0.008 

0.819 + 0.005 

0.890 + 0.007 

0.059 

0 

48 

NGVS-J123327.85 + 070649.5 

188.3660383 

7.1137599 

0.023 

20.978 + 0.002 

1.143 + 0.008 

0.752 + 0.005 

0.896 + 0.008 

0.168 

0 

49 

NGVS-J123359.09 + 074244.6 

188.4962125 

7.7124026 

0.022 

20.319 + 0.002 

1.124 + 0.005 

0.829 + 0.003 

0.969 + 0.004 

0.021 

1 

50 

NGVS-J123402.20 + 064950.3 

188.5091560 

6.8306518 

0.022 

21.218 + 0.003 

0.939 + 0.009 

0.651 + 0.006 

0.709 + 0.010 

0.084 

1 


NOTES - (1) Object ID number; (2) Object name in NGVS; (3) Right Ascension; (4) Declination; (5) Galactic extinction 
(Schlegel et ah 1998); (6) Aperture-corrected g magnitude within a 3-arcsec diameter aperture. Only SExtractor errors are 
accounted for here; (7-9) Colors measured in a 1.5-arcsec diameter aperture; (10) Envelope parameter defined in 8; (11) Object 
flag, 0 = not UCD, 1 = confirmed or possible UCD. 
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TABLE 5 

Structural Properties and Redshift Information for UCD Candidates in the M49 Region 


ID 

NGVSID 

9k 

'^k 

^h,g 

'^h,i 

t;^(MMT) 

Class 

u* gz 

UCD 



(mag) 

(mag) 

(pc) 

(pc) 

(km s“l) 




(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

1 

NG VS-J122708.94 + 074228.3 

24.228 + 0.010 

23.403 + 0.010 

45.60 + 0.70 

44.22 + 0.65 


1 

1 

1 

2 

NG VS-J122718.40 + 081154.2 

25.947+0.026 

25.062 + 0.047 

13.48 + 0.82 

18.41 + 1.46 


3 

1 

0 

3 

NGVS-J122735.11 + 072518.7 

26.010 + 0.025 

25.270 + 0.014 

15.53 + 0.89 

18.11 + 0.64 

38059 

1 

1 

0 

4 

NGVS-J122741.97+075613.3 

26.139 + 0.019 

25.337+0.024 

23.11 + 0.74 

23.14 + 0.85 

31816 

1 

1 

0 

5 

NGVS-J122744.26 + 072542.5 

25.688 + 0.006 

24.948 + 0.007 

12.77 + 0.37 

14.70 + 0.40 

1107 

1 

1 

1 

6 

NGVS-J122803.46 + 065343.1 

26.373 + 0.014 

25.524 + 0.037 

21.72 + 0.70 

28.08+1.46 


1 

1 

1 

7 

NGVS-J122806.14 + 065909.0 

25.353 + 0.014 

24.326 + 0.015 

15.98+1.12 

20.12 + 0.58 


1 

1 

1 

8 

NG VS-J122806.31 + 080342.8 

26.529 + 0.007 

25.861 + 0.012 

13.06 + 0.36 

7.84 + 0.48 

15078 

1 

1 

0 

9 

NG VS-J122818.37+074924.9 

26.586 + 0.014 

25.692 + 0.019 

21.60 + 0.58 

28.37 + 0.89 

41827 

4 

1 

0 

10 

NG VS-J122828.41 + 083115.5 

25.980 + 0.013 

25.221 + 0.007 

14.70 + 0.60 

19.45 + 0.49 

31226 

3 

1 

0 

11 

NG VS-J122828.60 + 070228.2 

25.771 + 0.008 

24.799 + 0.019 

20.71 + 0.49 

27.63+1.38 


1 

1 

1 

12 

NG VS-J122830.49 + 071929.1 

25.495 + 0.007 

24.763 + 0.016 

17.39 + 0.48 

12.68 + 0.94 

1528 

2 

1 

0 

13 

NG VS-J122849.25 + 075919.4 

26.269 + 0.006 

25.541 + 0.009 

14.34 + 0.30 

8.11 + 0.45 


1 

1 

1 

14 

NGVS-J122850.66 + 082621.0 

26.601 + 0.015 

25.829 + 0.010 

11.53 + 0.74 

14.39 + 0.31 

65197 

1 

1 

0 

15 

NGVS-J122857.99 + 083513.8 

23.721 + 0.020 

22.670 + 0.009 

28.03 + 0.91 

28.08 + 0.52 


3 

1 

0 

16 

NGVS-J122858.00 + 071759.9 

25.228 + 0.011 

24.205 + 0.017 

20.64 + 0.57 

21.98 + 0.68 

25640 

1 

1 

0 

17 

NGVS-J122906.90 + 073558.8 

26.222 + 0.018 

25.375 + 0.030 

16.55 + 0.80 

21.17+1.19 

51723 

1 

1 

0 

18 

NGVS-J122917.85 + 082558.1 

23.435 + 0.012 

22.505 + 0.022 

44.34 + 0.70 

36.59+1.34 

26908 

1 

1 

0 

19 

NG VS-J122922.85 + 074901.2 

25.782 + 0.010 

25.037+0.014 

15.13 + 0.39 

14.98 + 0.43 

1675 

1 

1 

1 

20 

NG VS-J122924.89 + 072105.3 

26.332 + 0.021 

25.666 + 0.020 

22.05 + 1.01 

19.92 + 0.98 


3 

1 

0 

21 

NGVS-J122925.10 + 075407.0 

25.185 + 0.004 

24.423 + 0.006 

15.91 + 0.24 

16.87 + 0.24 

1107 

1 

1 

1 

22 

NG VS-J122926.23 + 081658.8 

26.503 + 0.006 

25.739 + 0.010 

12.56 + 0.19 

8.23 + 0.31 

465 

1 

1 

1 

23 

NG VS-J122928.21 + 075355.4 

24.867+0.004 

23.960 + 0.005 

18.93 + 0.22 

16.62 + 0.24 

637 

1 

1 

1 

24 

NGVS-J122928.62 + 075824.6 

25.905 + 0.017 

25.140 + 0.018 

11.81 + 0.89 

14.55 + 0.51 

1114 

1 

1 

1 

25 

NGVS-J122934.41 + 080540.9 

25.478 + 0.021 

24.784 + 0.013 

20.92 + 0.97 

20.89 + 0.62 

1014 

4 

1 

0 

26 

NGVS-J122948.37+080042.0 

25.738 + 0.009 

24.931 + 0.013 

14.37 + 0.58 

15.33 + 0.58 

907 

1 

1 

1 

27 

NGVS-J122957.08 + 074820.0 

24.662 + 0.008 

23.916 + 0.013 

19.60 + 0.34 

14.55 + 0.34 

843 

1 

1 

1 

28 

NGVS-J122958.06 + 080409.2 

24.694 + 0.006 

23.787+0.010 

16.93 + 0.24 

18.49 + 0.36 

606 

1 

1 

1 

29 

NGVS-J123002.49 + 075846.7 

25.946 + 0.008 

25.211 + 0.011 

18.06 + 0.39 

17.50 + 0.43 

1820 

1 

1 

1 

30 

NGVS-J123004.35 + 073932.2 

26.586 + 0.007 

25.902 + 0.013 

12.20 + 0.48 

6.89 + 0.70 

1043 

1 

1 

1 

31 

NGVS-J123006.03 + 074901.2 

26.253 + 0.005 

25.543 + 0.008 

13.32 + 0.25 

13.09 + 0.40 

549 

1 

1 

1 

32 

NGVS-J123006.89 + 080509.5 

25.431 + 0.009 

24.674 + 0.011 

13.39 + 0.27 

13.47 + 0.30 

1001 

1 

1 

1 

33 

NGVS-J123008.39 + 082507.4 

24.074 + 0.010 

23.232 + 0.010 

17.47 + 0.58 

17.23 + 0.39 


1 

1 

1 

34 

NGVS-J123008.69 + 065353.4 

25.172 + 0.016 

24.223 + 0.011 

15.74+1.04 

25.41 + 0.58 

52650 

1 

1 

0 

35 

NGVS-J123009.17+074127.5 

24.567+0.006 

23.643 + 0.014 

12.25 + 0.34 

10.46 + 0.61 


1 

1 

1 

36 

NGVS-J123015.56 + 083445.0 

25.865 + 0.023 

24.954 + 0.048 

10.86 + 0.98 

11.37+1.44 

-129 

1 

1 

1 

37 

NGVS-J123016.65 + 080454.4 

25.815 + 0.006 

25.017+0.011 

16.93 + 0.34 

17.94 + 0.43 

778 

1 

1 

1 

38 

NGVS-J123018.06 + 075024.9 

24.352 + 0.004 

23.486 + 0.004 

31.19 + 0.22 

31.86 + 0.37 

1070 

1 

1 

1 

39 

NGVS-J123029.94 + 074755.5 

25.394 + 0.004 

24.567+0.006 

11.08 + 0.16 

11.16 + 0.22 

1395 

1 

1 

1 

40 

NGVS-J123047.31 + 073619.4 

25.158 + 0.006 

24.453 + 0.012 

14.60 + 0.33 

12.29 + 0.68 

1522 

2 

1 

0 

41 

NG VS-J123048.49 + 065420.2 

26.449 + 0.006 

25.611 + 0.023 

15.04 + 0.28 

18.48 + 0.71 

52822 

1 

1 

0 

42 

NGVS-J123109.91 + 080405.3 

24.819 + 0.009 

24.091 + 0.011 

14.79 + 0.34 

14.70 + 0.37 

426 

1 

1 

1 

43 

NG VS-J123113.76 + 074811.9 

26.239 + 0.008 

25.496 + 0.011 

17.26 + 0.48 

17.91 + 0.45 

1087 

1 

1 

1 

44 

NG VS-J123148.91 + 071620.9 

24.034 + 0.006 

23.124 + 0.014 

40.49 + 0.39 

46.69 + 0.98 


3 

1 

0 

45 

NGVS-J123152.37+064659.4 

25.003 + 0.015 

24.064 + 0.020 

36.99 + 0.82 

46.39+1.18 

46295 

1 

1 

0 

46 

NGVS-J123212.59 + 074047.6 

26.126 + 0.014 

25.310 + 0.023 

24.95 + 0.79 

31.38+1.04 

39582 

1 

1 

0 

47 

NGVS-J123320.46 + 080226.6 

25.856 + 0.018 

25.025 + 0.013 

27.26+1.01 

29.00 + 0.82 

36104 

1 

1 

0 

48 

NGVS-J123327.85 + 070649.5 

25.955 + 0.027 

25.129 + 0.035 

19.86+1.16 

23.93+1.40 


3 

1 

0 

49 

NGVS-J123359.09 + 074244.6 

25.292 + 0.008 

24.486 + 0.014 

13.96 + 0.39 

11.37 + 0.36 

854 

1 

1 

1 

50 

NGVS-J123402.20 + 064950.3 

26.175 + 0.009 

25.359 + 0.012 

25.52 + 0.51 

35.03 + 0.67 


1 

1 

1 


NOTE- (1) Object ID number; (2) Object name in NGVS; (3-4) King model magnitudes in g and i bands; (5-6) Half-light 
radii in g and i bands; (7) Radial velocity from MMT/Hectospec; (8) Class parameter: 1 = probable UCD, 2 = dwarf nucleus, 3 
= background galaxy, 4 = other type of contaminant; (9) Flag for whether the object selected in the u* gz color-color diagram: 
0 = no, 1 = yes; (10) Object flag, 0 = not UCD, 1 = conhrmed or possible UCD. 
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TABLE 6 

Photometric Properties of UCD Candidates in the M60 Region 


ID 

NGVSID 

«J2000 

<542000 

E{B - V) 

916 

(•u-* - 9)8 

(9 - i)s 

(9 - 

^env 

UGD 



(deg) 

(deg) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

1 

NGVS-J123909.32 + 112147.9 

189.7888329 

11.3633175 

0.033 

21.159 + 0.003 

1.415 + 0.010 

1.004 + 0.005 

1.175 + 0.007 

0.057 

1 

2 

NGVS-J123948.33 + 105157.0 

189.9513674 

10.8658447 

0.022 

20.858 + 0.002 

1.244 + 0.009 

0.907 + 0.004 

1.043 + 0.006 

0.130 

0 

3 

NG VS-J 123949.21 + 110135.3 

189.9550398 

11.0264758 

0.024 

21.167 + 0.003 

1.338 + 0.011 

0.981 + 0.005 

1.173 + 0.007 

0.078 

1 

4 

NGVS-J123950.99 + 112805.1 

189.9624727 

11.4680886 

0.039 

20.763 + 0.002 

1.208 + 0.007 

0.934 + 0.004 

1.162 + 0.005 

0.114 

0 

5 

NG VS-J124002.08 + 105517.2 

190.0086686 

10.9214418 

0.022 

21.420 + 0.003 

0.994 + 0.012 

0.664 + 0.007 

0.735 + 0.012 

0.261 

1 

6 

NG VS-J124049.54+121138.1 

190.2064072 

12.1939068 

0.039 

21.348 + 0.003 

1.219 + 0.011 

0.981 + 0.005 

1.127+0.008 

0.083 

1 

7 

NGVS-J124115.49 + 103811.9 

190.3145415 

10.6366475 

0.024 

21.350 + 0.003 

1.119 + 0.011 

0.865 + 0.006 

1.037+0.009 

0.069 

1 

8 

NGVS-J124151.11 + 105311.7 

190.4629499 

10.8865742 

0.021 

19.805 + 0.001 

1.405 + 0.005 

1.014 + 0.002 

1.242 + 0.003 

0.306 

0 

9 

NGVS-J124205.30 + 103959.1 

190.5220991 

10.6664070 

0.026 

21.462 + 0.003 

0.977+0.010 

0.729 + 0.007 

0.828 + 0.010 

0.070 

1 

10 

NGVS-J124206.89+115144.7 

190.5287016 

11.8624190 

0.028 

20.966 + 0.002 

1.113 + 0.007 

0.828 + 0.004 

0.963 + 0.006 

0.061 

1 

11 

NGVS-J124209.61 + 113654.7 

190.5400246 

11.6151888 

0.030 

20.940 + 0.002 

1.051 + 0.007 

0.782 + 0.004 

0.881 + 0.006 

0.059 

1 

12 

NGVS-J124216.65 + 114428.6 

190.5693756 

11.7412681 

0.031 

21.387 + 0.007 

1.468 + 0.023 

0.998 + 0.012 

1.258 + 0.012 

0.043 

1 

13 

NG VS-J124229.83 + 112029.7 

190.6242758 

11.3415871 

0.027 

19.975 + 0.001 

1.169 + 0.004 

0.830 + 0.002 

0.977+0.003 

0.045 

1 

14 

NG VS-J124232.60 + 115702.6 

190.6358256 

11.9507163 

0.028 

21.166 + 0.003 

1.012 + 0.008 

0.780 + 0.005 

0.896 + 0.008 

0.040 

1 

15 

NG VS-J124235.74 + 114254.2 

190.6489118 

11.7150453 

0.030 

21.423 + 0.003 

1.169 + 0.011 

0.888 + 0.006 

1.080 + 0.008 

0.055 

1 

16 

NGVS-J124238.75 + 115646.8 

190.6614439 

11.9463314 

0.028 

21.602 + 0.004 

0.910 + 0.011 

0.715 + 0.008 

0.786 + 0.012 

0.056 

0 

17 

NGVS-J 124244.72 + 111240.5 

190.6863239 

11.2112384 

0.027 

20.593 + 0.002 

1.246 + 0.007 

0.832 + 0.004 

1.007+0.005 

0.115 

1 

18 

NGVS-J124248.01 + 115513.8 

190.7000359 

11.9204879 

0.027 

19.812 + 0.001 

1.025 + 0.003 

0.818 + 0.002 

0.930 + 0.003 

0.120 

1 

19 

NGVS-J124312.47+111230.2 

190.8019424 

11.2083931 

0.030 

21.313 + 0.003 

1.148 + 0.011 

0.749 + 0.006 

0.888 + 0.009 

0.061 

1 

20 

NGVS-J124315.49+113922.8 

190.8145574 

11.6563462 

0.027 

20.719 + 0.002 

1.014 + 0.005 

0.666 + 0.004 

0.792 + 0.005 

0.037 

1 

21 

NGVS-J124324.30 + 112343.2 

190.8512567 

11.3953461 

0.024 

20.782 + 0.002 

0.987+0.006 

0.732 + 0.004 

0.814 + 0.006 

0.083 

1 

22 

NGVS-J124352.42 + 112534.2 

190.9684091 

11.4261669 

0.024 

18.830 + 0.001 

1.412 + 0.002 

0.968 + 0.001 

1.165 + 0.001 

0.080 

1 

23 

NGVS-J 124352.48 + 112518.3 

190.9686634 

11.4217508 

0.025 

20.312 + 0.001 

1.072 + 0.005 

0.773 + 0.003 

0.875 + 0.004 

0.048 

1 

24 

NGVS-J 124437.85 + 103801.1 

191.1577272 

10.6336481 

0.026 

20.012 + 0.001 

1.110 + 0.004 

0.858 + 0.002 

0.989 + 0.003 

0.085 

1 

25 

NGVS-J 124524.02 + 110535.2 

191.3500929 

11.0930994 

0.031 

19.454 + 0.001 

1.397+0.003 

1.034 + 0.002 

1.158 + 0.002 

0.055 

1 

26 

NGVS-J 124531.21 + 112357.2 

191.3800543 

11.3992236 

0.027 

20.351 + 0.001 

1.083 + 0.005 

0.806 + 0.003 

0.898 + 0.004 

0.112 

1 

27 

NGVS-J 124532.99 + 114125.3 

191.3874620 

11.6903730 

0.026 

21.275 + 0.003 

1.011 + 0.009 

0.809 + 0.006 

0.978 + 0.008 

0.138 

0 

28 

NGVS-J 124538.98 + 114526.8 

191.4124258 

11.7574345 

0.028 

21.282 + 0.004 

1.034 + 0.013 

0.702 + 0.009 

0.786 + 0.013 

0.118 

0 

29 

NGVS-J124539.33+104808.1 

191.4138693 

10.8022486 

0.024 

19.993 + 0.001 

1.239 + 0.005 

0.899 + 0.003 

1.041 + 0.003 

0.099 

1 


NOTES - (1) Object ID number; (2) Object name in NGVS; (3) Right Ascension; (4) Declination; (5) Galactic extinction 
(Schlegel et ah 1998); (6) Aperture-corrected g magnitude within a 3-arcsec diameter aperture. Only SExtractor errors are 
accounted for here; (7-9) Colors measured in a 1.5-arcsec diameter aperture; (10) Envelope parameter defined in 8; (11) Object 
flag, 0 = not UCD, 1 = confirmed or possible UCD. 


TABLE 7 

Structural Properties and Redshift Information for UCD Candidates in the M60 Region 


ID 

NGVSID 

9k 


'^h,g 


T'r(AAT) 

Glass 

u* gz 

UCD 



(mag) 

(mag) 

(pc) 

(pc) 

(km s“^) 




(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 


1 

NGVS-J123909.32+112147.9 

26.190 + 0.009 

25.224 + 0.011 

12.39 + 0.68 

12.60 + 0.58 


1 

1 

1 

2 

NGVS-J123948.33+105157.0 

25.805 + 0.029 

24.874 + 0.039 

21.08 + 0.94 

23.58+1.43 


3 

1 

0 

3 

NGVS-J123949.21 + 110135.3 

26.214 + 0.009 

25.268 + 0.009 

11.72 + 0.43 

11.56 + 0.30 


1 

1 

1 

4 

NGVS-J123950.99+112805.1 

25.740 + 0.024 

24.742 + 0.012 

21.20+1.06 

27.76 + 0.80 


3 

1 

0 

5 

NGVS-J124002.08 + 105517.2 

26.523 + 0.020 

25.875 + 0.025 

12.72 + 0.57 

14.98 + 0.77 


1 

1 

1 

6 

NGVS-J124049.54+121138.1 

26.376 + 0.008 

25.463 + 0.021 

13.93 + 0.45 

13.44+1.13 


1 

1 

1 

7 

NGVS-J124115.49+103811.9 

26.359 + 0.013 

25.518 + 0.008 

14.31 + 0.64 

13.39 + 0.51 


1 

1 

1 

8 

NGVS-J124151.11 + 105311.7 

24.887 + 0.017 

23.810 + 0.013 

20.24 + 0.82 

25.90 + 0.65 


3 

1 

0 

9 

NGVS-J124205.30+103959.1 

26.438 + 0.020 

25.639 + 0.022 

18.14 + 0.71 

22.54 + 0.83 


1 

1 

1 

10 

NGVS-J124206.89+115144.7 

25.969 + 0.008 

25.164 + 0.010 

11.87+0.27 

10.53 + 0.43 


1 

1 

1 

11 

NGVS-J124209.61 + 113654.7 

25.920 + 0.004 

25.147 + 0.010 

12.97+0.18 

11.69 + 0.33 


1 

1 

1 

12 

NGVS-J124216.65 + 114428.6 

26.407 + 0.012 

25.425 + 0.013 

11.56 + 0.51 

10.47+0.79 


1 

1 

1 

13 

NGVS-J124229.83+112029.7 

24.979 + 0.010 

24.162 + 0.009 

14.25 + 0.31 

15.79 + 0.36 

1493 

1 

1 

1 

14 

NGVS-J124232.60 + 115702.6 

26.158 + 0.005 

25.389 + 0.007 

12.97+0.22 

12.69 + 0.42 


1 

1 

1 

15 

NGVS-J124235.74+114254.2 

26.453 + 0.010 

25.573 + 0.014 

11.34 + 0.58 

10.53 + 0.85 


1 

1 

1 

16 

NGVS-J124238.75+115646.8 

26.600 + 0.029 

25.805 + 0.036 

16.95 + 0.98 

21.98+1.18 


3 

1 

0 

17 

NGVS-J124244.72+111240.5 

25.462 + 0.029 

24.536 + 0.016 

31.07+1.18 

38.77+1.28 


1 

1 

1 

18 

NGVS-J124248.01 +115513.8 

24.847 + 0.007 

24.048 + 0.008 

13.94 + 0.28 

12.60 + 0.61 

1241 

1 

1 

1 

19 

NGVS-J124312.47+111230.2 

26.323 + 0.009 

25.490 + 0.010 

15.59 + 0.51 

20.80 + 0.48 


1 

1 

1 

20 

NGVS-J124315.49+113922.8 

25.697 + 0.011 

25.020 + 0.018 

11.43 + 0.42 

11.64 + 0.51 


1 

1 

1 

21 

NGVS-J124324.30+112343.2 

25.684 + 0.007 

24.986 + 0.010 

24.82 + 0.40 

24.98 + 0.52 


1 

1 

1 

22 

NGVS-J124352.42+112534.2 

23.872 + 0.008 

22.911 + 0.008 

13.94 + 0.36 

14.66 + 0.37 


1 

1 

1 

23 

NGVS-J124352.48 + 112518.3 

25.249 + 0.008 

24.496 + 0.011 

17.51 + 0.42 

17.75 + 0.58 


1 

1 

1 

24 

NGVS-J124437.85+103801.1 

24.953 + 0.007 

24.139 + 0.007 

19.83 + 0.30 

18.41 + 0.28 


1 

1 

1 

25 

NGVS-J124524.02 + 110535.2 

24.481 + 0.009 

23.521 + 0.018 

11.26 + 0.42 

10.85 + 0.51 


1 

1 

1 

26 

NGVS-J124531.21 + 112357.2 

25.419 + 0.011 

24.630 + 0.010 

11.10 + 0.40 

11.17+0.52 


1 

1 

1 

27 

NGVS-J124532.99+114125.3 

26.300 + 0.022 

25.351 + 0.028 

17.45+1.00 

26.01 + 1.32 


3 

1 

0 

28 

NGVS-J124538.98+114526.8 

26.192 + 0.030 

25.446 + 0.013 

25.03+1.22 

27.18 + 0.54 


3 

1 

0 

29 

NGVS-J124539.33+104808.1 

25.024 + 0.020 

24.154 + 0.009 

11.53 + 0.76 

16.07+0.49 


1 

1 

1 


NOTE- (1) Object ID number; (2) Object name in NGVS; (3-4) King model magnitudes in g and i bands; (5-6) Half-light 
radii in g and i bands; (7) Radial velocity from 2dE/AAOmega AAT program (Zhang et al. 2015); (8) Class parameter: 1 = 
probable UCD, 2 = dwarf nucleus, 3 = background galaxy, 4 = other type of contaminant; (9) Elag for whether the object 
selected in the u*gz color-color diagram: 0 = no, 1 = yes; (10) Object flag, 0 = not UCD, 1 = confirmed or possible UCD. 
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TABLE 8 

UCD Samples in Program Regions 


Region 

(1) 

(aj 20 oo) 

(deg) 

(2) 

(<5j200o) 

(deg) 

(3) 

n 

(deg^) 

(4) 

MjCD.ugz 

(5) 

SuCD.ugz 

(deg^) 

(6) 

■^UCD.uiKs 

(7) 

SuCD.uiKg 

beg^) 

(8) 

M87 

187.60 

12.44 

3.61872 

no 

30.4±2.9 

96 

26.5±2.7 

NGVS+0+0 

187.94 

12.05 

0.90390 

50 

55.3±7.8 

47 

52.0±7.6 

NGVS+0+1 

188.03 

12.94 

0.90314 

17 

18.8±4.6 

17 

18.8±4.6 

NGVS-1+0 

187.15 

12.01 

0.90656 

14 

15.4±4.1 

6 

6.6±2.7 

NGVS-1+1 

187.09 

12.94 

0.90512 

29 

32.0±5.9 

26 

28.7±5.6 

M49 

187.57 

07.75 

3.61230 

40 

11.1±1.8 



NGVS-1-4 

187.14 

08.20 

0.90095 

19 

21.1±4.8 



NGVS-1-5 

187.11 

07.28 

0.90525 

11 

12.2±3.7 



NGVS+0-4 

188.04 

08.20 

0.90000 

5 

5.6±2.5 



NGVS+0-5 

188.05 

07.27 

0.90610 

5 

5.5±2.5 



M60 

190.50 

11.53 

3.61201 

23 

6.4±1.3 



NGVS+3+0 

190.93 

11.98 

0.90546 

7 

7.7±2.9 



NGVS+3-1 

190.93 

11.10 

0.90080 

11 

12.2±3.7 



NGVS+2+0 

189.99 

12.01 

0.90265 

1 

l.lil.l 



NGVS+2-1 

189.98 

11.03 

0.90310 

4 

4.4±2.2 



Background Eields 



4.04765 

9 

2.2±0.7 



NGVS_BG1 

202.10 

21.03 

1.01687 

2 

2.0±1.4 



NGVS_BG2 

181.12 

27.39 

1.01351 

1 

l.Oil.O 



NGVS_BG3 

171.75 

15.64 

1.00172 

4 

4.0±2.0 



NGVS_BG4 

203.39 

06.76 

1.01555 

2 

2.0±1.4 



GEHTLS-D2+UltraVISTA 

150.12 

02.24 

1.05779 



1 

0.9±0.9 


NOTE - (1) Region designation (2-3) Mean Right Ascension and Declination of the region; (4) Eield of view; (5-6) Number 
and number density of LCDs selected based on the gz color-color diagram; (7-8) Number and number density of LCDs 
selected based on the u^iKs color-color diagram. 








TABLE 9 

UCD Color Distributions in the M87 region 
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